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Liquid=Solid Equilibrium in Copper and Silver Alloys; with an 
Appendix on the Equilibrium Diagram of the System 
Copper—Arsenic} 


By Witu1am Hume-Roruery and JENNIFER BuRNS 
Department of Metallurgy, Oxford 


[Received June 17, 1957] 


SUMMARY 

The liquidus and solidus curves for the «-solid solutions, of Cd, In, Sn, 
and Sb in Ag, and of Zn, Ga, Ge, and As, in Cu are discussed. In the 
Ag-series a given depression of the solidus is produced by atomic per- 
centages of Cd, In, Sn, and Sb proportional to 1/2? : 1/32: 1/4? : 1/62, 
and in the Cu-series by amounts of Zn and Ga proportional to 1/3? : 1/42. 
The increases in mean lattice spacing produced by a given atomic per- 
centage of Cd, In, Sn, and Sb dissolved in Ag are as 2:3: 4:6, and by a 
given atomic percentage of Zn and Ga dissolved in Cu are as 3:4. The 
solidus curves and lattice distortions thus involve the same whole number 
sequences, but this correspondence breaks down in the system Cu—As 
where the electrochemical factor is high, and the same effect to a lesser 
degree is apparent in the system Cu-Ge. In the systems Ag—Cd, Ag—In, 
Ag-Sn, Cu-Zn, Cu-Ga, and Cu—Ge the liquidus and solidus curves in 
dilute solutions are connected by a relation of the type xj=k,/x, where 
x, and x, are the atomic percentages of solute in the liquid and solid 
phases. In this class of alloy the liquidus and solidus curves depend 
primarily upon electron concentration, modified by lattice distortion, 
and diagrams are given showing these effects. Increasing electrochemical 
factor (e.g. Cu-As, Ag—Mg) produces an additional steeper fall in the 
liquidus and solidus. The equations of H. Jones are used to calculate 
the effects of solute elements on the latent heat of fusion of Ag and Cu, 
and the way in which these are affected by electron concentration and 
lattice distortion is discussed. 


§ 1. INTRODUCTION 


Tue solid solutionst of Cd, In, Sn, and Sb in silver, and of Zn, Ga, Ge, 
and As in copper show many properties for which the changes produced 
by a given atomic percentage of solute increase with the valency (v) 
of the latter. The increases in the mean lattice spacing produced by a 
given atomic percentage of Cd, In, Sn, and Sb dissolved in silver are as 


+ Communicated by the Author. ik ES, 
+ Chemical symbols are used for abbreviation, and for clarity in describing 
the relative effects of different elements. 
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2:3:4:6 (Hume-Rothery and Reynolds 1937, Owen and Roberts 
1939), whilst for Zn, Ga, Ge, and As dissolved in copper the increases are 
as 3:4:5:7. In the liquid=solid equilibria for the «-solid solutions in 
silver, a given depression of the liquidus is produced by amounts of Cd, 
In, Sn, and Sb proportional to 1/2a:1/3q@: 1/46: 1/56, where a and 6 are 
very nearly equal (Hume-Rothery, Mabbott, and Channel-Evans 1934f 
Hume-Rothery and Reynolds 1937 a), and H-R, M, and C-E showed that 
to a first approximation a similar relation existed for the «-solid solutions 
of Zn, Ga, and Ge in copper. These authors also claimed that the 
corresponding depressions of the solidus curves were roughly propor- 
tional to 1/v?, but this relation was not exact. 

The curious whole number relations for the lattice distortions and 
depressions of liquidus curves are not understood, and the experimental 
accuracy is not sufficiently great to preclude a departure of 0-1 from the 
exact whole numbers. Since the liquidus and solidus curves give the 
compositions of the two phases in equilibrium with each other, any really 
fundamental relation will probably involve the compositions of both 
phases, and although a definite valency effect is clear, the apparent 
agreement with whole number sequences may be of no significance. 

The general theory of alloys suggests that the structures of this type 
of alloy are determined by the interplay of three factors (1) electron 
concentration, (2) size-factor, and (3) electrochemical factor, the last- 
named being important only when the solvent and solute differ appreci- 
ably in electrochemical characteristics. In the present paper we discuss 
the liquidus and solidus curves in greater detail, and give new results 
for the system copper—arsenic in which the electrochemical factor is 
considerable. The discussion is confined to alloys of copper and silver 
with solute elements whose valencies are higher, and whose ionic radii 
are smaller than those of the solvent. Under such conditions the lattice 
distortion is due mainly to the valency electrons of the solute. 


§ 2. THe EQuinisrium Diagrams 


| Figure | shows the equilibrium diagrams in such a way as to distinguish 
between the temperature-composition limits in which the alloys are 
totally solid or totally liquid. In the systems Ag—Cd and Cu—Zn, the 
liquidus curves show only slight deviations from a smooth curve hotween 
the melting points of the two metals, and except at the extreme Cd-rich 
and Zn-rich ends of the diagram the horizontal portions of the solidus 
curves are short. In the systems Ag—In and Cu-Ga, the liquidus curves 
are again nearly smooth curves across the diagram, but the solidus curves 
show marked falls after about 35 atm°% solute. 

The system Ag—Sn shows the same characteristics except that the 
marked widening in the freezing range occurs at about 25 atm% Sn 
In contrast to this, the system Cu-Ge shows a slight canna = 


+ This paper will be referred to as H-R, M, and C-E, 
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26 atm% Ge, and the generally falling liquidus and solidus curves end in 
a eutectic at 37 atm% Ge. 

On passing to the elements of Group V, this difference is accentuated. 
In the system Ag-Sb there is a general fall in the liquidus and solidus 
until a eutectic is reached in the middle of the diagram. In contrast 
to this, the system Cu—As shows a pronounced maximum at 26 atm% As. 


Fig. 1 


In these diagrams the upper lines are the liquidus curves of the systems con- 
cerned. The shaded areas are the (liquid-+solid) 2-phase fields, and 
the areas below these are those in which the alloys are totally solid. 


We may interpet these differences as showing that, in the systems 
Ag—Cd, Ag—In, Ag—Sn, Cu-Zn, and Cu-Ga, the electrochemical factor is 
comparatively small, so that the equilibrium diagrams are determined 
mainly by the electron concentration and size-factor (i.e. lattice distor- 
tion) effects. In the system Cu-Ge, the electrochemical factor is already 
beginning to affect the form of the diagram, and is playing a far more 
important part in the system Cu—As. This interpretation is confirmed 
by the fact that the equilibrium diagram of the system Cu-Si is very 
like that of the system Cu—Ge, but with a much more pronounced maxi- 
mum in the liquidus curve, as would be expected from the more electro- 


negative nature of silicon. 


§ 3. Tue Liquipus AND Souipus Curves: Empirical RELATIONS 

In view of recent data (see p. 1195) H-R, M, and C-K’s claim that the 
depression of the solidus curves is proportional to 1/v® is incorrect. 
Instead of this we find that for dilute solutions in the silver series, a 
given depression of the solidus curves is produced by atomic percentages 
of Cd, In, Sn, and Sb proportional to 1/2? : 1/3? : 1/4: 1/6?, whilst in 
the copper series a given depression of the solidus curves is produced 
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by amounts of Zn, and Ga proportional to 1/32: 1/42. These whole 
numbers are the same as those involved for the lattice distortionsf. 

We also find that although H-R, M, and C-E’s claim that a given 
depression of the liquidus curves was produced by atomic percentages of 
Zn, Ga, and Ge proportional to 1/2: 1/3: 1/4 was a rough general 
approximation, the relative depressions for Zn and Ga are in much better 
agreement with the factors 1/3 : 1/4, where the whole numbers are again 
those involved for the lattice distortions. Even though the whole 
number relations may not be exact, the general correspondence with 
the mean lattice distortions is undoubted. 

These relations imply that, for the systems Ag-Cd, Ag—In, Ag—Sn, 
Cu-Zn, and Cu-Ga, at a given temperature the concentrations 2 and 
a, of the solute in the liquid and solid phases are connected by a relation of 
the type 1=k,/x,. This is the same relation as that given by H-R, 
M, and C-E, although the whole numbers suggested by these authors 
for Cu-Zn and Cu-Ga were wrong. Such a relation will not hold for 
systems in which the solid solution is very restricted, because the liquidus 
and solidus are connected thermodynamically so that 2 retains a definite 
value when 2,=0, ie. when no solid solution is formed. For the case of 
monatomic dilute solutions the value of x for 7,0 can be calculated 
by means of Van’t Hoff’s equation 


Al 1 
an an pm): 

In fig. 2, x for the silver alloys at 900°c is plotted against 1/2,, together 
with the value for x,=0 calculated from the above equation. The 
points for Ag—Cd, Ag—In, and Ag-Sn lie on a straight line through the 
origin. The points for Ag—-Zn and Ag—Hg, where the divalent solutes 
are from adjacent Periods also lie on this straight line, although at 
900°c the solidus curve in the system Ag—Hg involves such a long extra- 
polation from the observed points that the test is not a critical one. 
The points for Ag-Sb, Ag-Ga, Ag—Ge, and Ag—Pbt where a, is small show 
a gradual tailing off towards the calculated value for x,=0. The point 
for Ag—Al is in fair agreement with the straight line, and the only wide 
deviation is for the system Ag—Mg, where the electrochemical factor is 
large, and the equilibrium diagram shows a pronounced maximum in the 
iene region of the 8 phase, and the point lies above the line in 

(are 

The general conclusion to be drawn from fig. 2 is that for this class 
of alloy the curve connecting x, and x, approaches asymptotically to the 


calculated value for 2=0 in cases where x, is small, and to a straight 
eee ee ee eee 

} As explained before, the experimental accuracy does not entitle us to 
claim an exact whole number effect, but we use the whole numbers for abbrevia- 
tion in discussion. 

t In this system Ag-Pb the ionic radius of the Pb!* ion is smaller than that 
of Ag*, although if Pb dissolves in the incompletely ionized form the ionic 
radius of the Pb?* ion would be greater. 
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line through the origin (i.e. a direct proportionality) when zs is large 
provided that the electrochemical factors are small. 


2 


Fig. 2 


| : er 


In this diagram 2 is plotted against +/x,, where x, and 2, are the atomic 
percentages of solute in the liquid and solid («) phases in equilibrium 
with one another at 900°C in the alloys of silver with the elements 
shown. The value for z,=O0 is calculated by Van’t Hoff’s equation. 
Similar diagrams are obtained at any constant temperature above 900°c 
whilst at lower temperatures the relation gradually becomes less exact. 


Fig. 3 


all 
| 2 3 Fee ES 


2, plotted against 1/x, for copper alloys at 1035°c. See fig. 2. 


In fig. 3, x is plotted against 1/x, for the copper alloys at 1035°c ; 
this temperature is chosen because it is just above that of the 
(a+S—=Liquid) eutectic in the system Cu—Al and is thus the lowest 
temperature for which a point from this system can be included. The 
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points for Cu-Zn, Cu-Ga, and Cu-Ge lie on a straight line through the 
origin, whilst that for Cu—As lies above this line, as for Ag—-Sb in fig. 2. 
In the system Cu—Al, the depression of liquidus and solidus is abnormally 
small, for reasons which are not understood, and the Cu—Al point is above 
the line in fig. 3, as is the case for the Ag—Mg point in fig. 2, and Al and 
Mg are the two most electropositive of the solutes. The points for 
Cu-In and Cu-Sn where the solid solubility is small approach the value 
calculated for x,=0, but for these the ionic radii of the solutes are of the 
same order as that of Cur. 


§ 4, ELECTRON CONCENTRATION AND LaTTIce DISTORTION 


We may accept the view that, when the electrochemical factor is small, 
the liquidus and solidus curves are determined first by the electron 
concentration, modified by the size factor or lattice distortion. The 
mean lattice distortion as measured by ordinary x-ray diffraction is 
(apart from lattice vacancies) directly related to the change in atomic 
volume, but it is not a very satisfactory quantity because it fails to take 
into account the local distortions around the solute atoms. Attempts. 
to measure the local distortions have been made for other alloys by 
Warren, Averbach and Roberts (1951) but the results are not yet 
sufficiently precise for the present purpose. In a qualitative way some 
of this difficulty can be overcome by comparing mean lattice distortions. 
in alloys containing the same concentration of solute atoms, since in this 
case the distortion results from the same number of centres. This policy 
is adopted on p. 1188 whilst in other cases we can only use the mean 
values at different concentrations. In the absence of high-temperature 
data we can only use lattice-spacings determined at room temperature. 
This seems justified in view of the fact that Owen and Roberts (1939) 
showed that the relative effects of different elements on the lattice- 
spacings of these alloys were the same at 300°c. as at 18°c. 

In the alloys under consideration a given electron concentration is 
produced by atomic proportions of (Zn, Cd), (Ga, In), (Ge, Sn), and 
(As, Sb) proportional to 12:6:4:3. In the silver series where the 
mean lattice distortion produced by equiatomic percentages of Cd, In, 
Sn and Sb are as 2:3: 4:6, the lattice distortions at a given electron 
concentration are in the ratios : 


(122) Ag Cd: (6X3) Ag In: (4x4) Ag Sn: (3x6) Ag Sb, i.e. 
24°18 5 16718; 

In the copper series where the lattice distortion factors are 3, 4, 5 and 7, 
the lattice distortions at a given electron concentration are in the ratios : 
(12 3) Cu Zn : (6x 4) Cu Ga : (4 5) Cu Ge: (37) Cu As, ie. 

36 224 2202 
In each series, therefore, the mean lattice distortion for a given electron. 
concentration decreases markedly on passing from the divalent to the 


Ce oe 
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trivalent solute, with a further decrease on passing to the tetravalent 
solute, after which there is a slight increase on passing to the solute of 
Group V. 

§ 5. Tue «-Liguipus anp Sourpus CurvEs 

Figure 4 shows the liquidus and solidus curves for the above series of 
alloys plotted in terms of electron concentration, and the following 
points may be noted—for abbreviation we shall refer to the 2-phase 
(Liquid+-Solid) region as the (L-++8) field. 

For the silver series, the width of the (L-+-S) fields increases with the 
valency of the solute. The (L+8) field for Ag—In is not placed sym- 
metrically about that for (Ag—Cd but is shifted clearly in the direction of 
higher temperature. On passing to Ag—Sn the (L-++8) field widens further, 
and is not quite symmetrical about the Ag—In (L+S8) field, and although 
slight, the shift is again in the direction of higher temperature. On 
passing to Ag-Sb this tendency is reversed and the greatly widened 
(L-+-S) region is moving in the direction of lower temperature. These 
changes run parallel to the factors for mean lattice distortions at con- 
stant electron concentrations which are (p. 1182) as 24:18:16: 18. 


Fig. 4 


1000°C 


900 


850 . = 
1.00 1.05 1.10 1.15 I 20! “E.G: 


«-liquidus and solidus curves plotted against electron concentration for Ag—Cd, 
Ag—In, Ag-Sn and Ag—Sb alloys. 


For the copper series (fig. 5) the relative widths and positions of the 
(L-LS) fields for the systems Cu-Zn and Cu—Ga resemble those for Ag—Cd 
and Ag—In, but for Cu-Ge and Cu—As this correspondence breaks down. 
Tn terms of electron concentration, the liquidus curves for Cu-Ga, Cu-Ge, 
and Cu-As are roughly similar, but the solidus curves for Cu-Ge and 
Cu—As are distinctly lower. The effect is very marked for Cu—As, and 
this suggests that increasing electrochemical factor acts so as to lower 
the general direction of the (L+-S) field. 

It was shown by Andrews and Hume-Rothery (1941) that in the 
a/8 brass kind of equilibrium a simple relation existed between the lattice 
distortions and the electron concentrations of the lines half way between 
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the compositions of the two phases in equilibrium. We may, therefore, 
define a quantity @ by the relation @=4(7,+-21) where x, and 2, are the 
atomic fractions of the solute in the solid and liquid phases at a given 


temperature. ; 
Fig. 5 


Meeyer 


1050 


1000 


950 | 
100 1-05 HO HS L20EC. 


«-liquidus and solidus curves plotted against electron concentration for Cu-Zn, 
Cu-Ga, Cu-Ge and Cu—As alloys. 


Fig. 6 
{000°C j= : 
Ag Cd and Ag Sn 


© O > OAgsb 
950 


900 


850 


1.0 Il 1.2 13 EG 4 
x versus electron concentration for Ag-Cd, Ag—In, Ag—-Sn, and Ag-Sb alloys. 
xis equal to 3(%,-+-a,) and is, thus, the line half-way between the liquidus 

and solidus curves. 


Figure 6 shows the relation between (in electron concentration) and tem- 
perature for the silver alloys. On passing from Ag—Cd > Ag—-In > Ag-Sn 
the @ curves rise, and then on passing to Ag-Sb the @ curve falls and 
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is almost the same as that for Ag—In. These effects are in the same 
order as the mean lattice distortions at a given electron concentra- 
tion, and the coincidence of the </Z7'’ curves for Ag—In and Ag—Sb is 
remarkable in view of the equality of the mean lattice distortions for a 
given electron concentration for these two solutes. 

Figure 7 shows the </7' relations for the copper alloys. For Cu-Zn 
and Cu-Ga the relative positions of the two curves resemble those for 
Ag—Cd and Ag—In in fig. 6, but on passing to Group IV the @/7' curve is 
lowered for Cu-—Ge, in contrast to the behaviour of the Ag-Sn curve in 
fig. 6, whilst the curve for Cu—As shows a striking fall, in agreement with 
the increasing electrochemical factor. 


Fig. 7 


1000 


900 
110 1:20 1:30 P40 EC. 


% versus electron concentration for Cu-Zn, Cu-Ga, Cu-Ge, and Cu—As alloys. 
See fig. 6. 


We may, therefore, conclude that the «/7' curves are determined prim- 
arily by electron concentration, and are depressed by increasing lattice 
distortion and by increasing electrochemical factor. We can thus 
understand why the #/7’ curve for the system Cu-Si nearly coincides 
with that for Cu-Ge. For in the Cu-Si «-solid solution, the lattice 
distortion is very small, and would favour a high #/7' curve, whilst the 
more electronegative nature of silicon would favour a low 2/7’ curve, 
and we can understand how the two effects tend to neutralize each 
other. 

§ 6. THE JU VaLuEs 

The liquidus and solidus curves discussed above are not independent, 
but are connected thermodynamically. Following Jones (1937) we 
may write F; and Fs for the free energies per atom for the liquid and 
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solid phases in equilibrium with one another at a given temperature, 
and let a and 2, be the compositions, i.e. the atomic fractions of the: 
solute in the two phases. The conditions for equilibrium are expressed. 
by the two equations : 

OF, OF,  F(%)—F (as) 


021. OX, «= LI —A 


(1), 


The free energy of each phase may be written in the form 
Fi,=U(@)+O+D. 2). tn 3 ep 


U(«) is the internal energy of the alloy at the absolute zero of tempera- 

ture as a function of its composition, measured with respect to the 

energies of the individual atoms. C is the free energy associated with the 

rise of temperature from the absolute zero to temperature 7’, and is thus. 

equal to: 
T gen 

{ CaT—T| zea. 
0 0 


D is the disorder term arising from the entropy of mixing and is equal to. 
kT {x In «+(1—z2) In (1—2)}. 


Jones assumed that the change in vibrational energy with composition 
at a fixed temperature was the same for both liquid and solid phases, 
so that the C term of eqn. (2) could be ignored in considering differences. 
between the free energies in liquid and solid states. In this way, Jones. 
obtained an expression for the quantity 4U defined by : 


AU=(U,—U)),—(U,—U1)y 


where (U,—U}); is the difference between the internal energies of the 
solid and (supercooled) liquid phases of composition % at the absolute 
zero, and (U,—U}), has the same meaning for the pure solvent. This. 
equation is fT : 
(lok las, Ga( 1m) 

au=( T ) A-key In =. +2 ne as (3): 
where \ and 7’, are the latent heat of fusion, and melting point of the 
pure solvent. 2, and x are the compositions of solid and liquid phases 
in equilibrium at temperature 7, and €=4(x,-++2) is the same quantity 
referred to on p. 1187. 

The second term in (3) arises from the disorder term in eqn. (1), and 
is small (5-10°%,) compared with the first term in (3). The increase in 
AU is thus determined mainly by the depression of the </7' curve. 

AU may be regarded as the decrease in the latent heat of fusion 
resulting from the addition of a fraction # of solute—more strictly it is 
the decrease in the latent heat of fusion of the solid to (hypothetical) 
supercooled liquid at the absolute zero. Positive values of JU mean 


+ This equation is printed wrongly in Jones’ paper, but the calculated 
values are correct. We have to thank Professor Jones for discussion of this. 
point. 
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that the liquidus and solidus are depressed by the addition of solute, 


whilst negative values of 4U correspond to a rise in melting point. 
The Jones equation, thus, enables the JU values to be calculated 


from the liquidus and solidus curves. Since the #/7' curves fall pro- 
gressively more steeply with increasing valency of the solute, the curves 


Fig. 8 

200 
TANG) 
CuZn CuGe CuGa 


100 


Agin AgSn 


ee te 


100 


° 
a “ul I2 EC. 


AU versus electron concentration for silver and copper alloys. For the 
definition of 4U, see p. 1186. 


connecting AU with # (in atomic concentrations) show, for a given value 
of € a progressive increase in 4U on passing along the series Cu-Zn, 
Ga, Ge and Ag—Cd, In, Sn, Sb, and Jones claimed that JU was very 
nearly proportional to the electron concentration. 
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Since the work of Jones more accurate data have become available, 
the details of which are given in Appendix B. Figure 8 shows the electron 
concentration plotted against the 4U values calculated from the new 
data, and the following points may be noted. 

(1) When plotted against electron concentration, the 4U curves are 
not superposed, and consequently 4U does not depend on electron con- 
centration alone. 

(2) In the copper series the 4U/e.c. curvest for Cu-Zn and Cu-Ga 
are displaced so that JU increases with the lattice distortion factors 
for a given electron concentration (p. 1185). For Cu—As, the 4U/e.c. curve 
is much higher than would be expected from the lattice distortion factors 
which as shown on p. 1182 are as 36: 24: 20: 21 for a given electron con- 
centration in the solutions of Zn, Ga, Ge, and As in copper. We may 
reasonably ascribe this to the electrochemical factor, and the fact that 
the curve for Cu Ge is above that for Cu Ga may be the result of the 
electrochemical factor beginning to exert its effect in the system Cu Ge 
(see p. 1179). 

(3) In the silver series where the lattice distortions at constant electron 
concentration are as 24:18: 16:18 for Ag—Cd, In, Sn, and Sb, the 4U 
curves for the system Ag—Cd, Ag—In, and Ag—Sn follow the lattice distor- 
tion sequence, but the Ag—Sb curve is higher than would be expected 
in contrast to the behaviour of the empirical # curves of fig. 6. This 
might suggest that in the system Ag—Sb, the electrochemical factor was 
beginning to exert its effect, but the approximations in Jones’ treatment 
do not justify further discussion. 


§ 7. THE RELATION BETWEEN 4U anp THE M@an LatrTicE 
DISTORTION 


As explained previously (p. 1182) the mean lattice distortion is directly 
related to the change in atomic volume, but fails to take into account 
the different local distortions which may occur round different solute 
elements. If we consider alloys containing a given atomic percentage of 
solute, the mean lattice distortions proceed from the same number of 
disturbing centres, and the conditions may be more directly comparable 
than if the atomic percentage of solute is a variable. 

Figure 9 shows AU plotted against the mean lattice distortions, Ja, 
for the Ag—Cd, In, Sn, Sb alloys containing 2, 6, and 12 atm°%, of solute. 
The points lie on straight lines, the extrapolations of which indicate 
negative values of JU for da=0. In all of these alloys the formation 
of the solid solution results in an increase in the electron concentration 
which is responsible for most of the lattice expansion. We may perhaps 
interpret the negative value of AU as indicating that, if in this series of 
alloys 4a were zero, there would have to be some strong attractive force 


{ For abbreviation we refer to the 4U versus electron concentration curves 
as the 4U/e.c. curves. 
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in order to neutralize the expansion resulting from the increased electron 
concentration, and this attraction would stabilize the solid phases and 
raise the melting point (iie—dU). In the system Ag-Sb, #@=12 
corresponds to a point slightly above the peritectic horizontal of the 
a-phase (fig. 1) and since the «-solid solution extends only to 7 atm Sb, 
we have assumed that the lattice distortion follows the 2:3:4:6 
sequence referred to on p. 1178. 


Fig. 9 
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AU is plotted against da for constant atomic percentages of Cd, In, Sn, and 


Sb dissolved in Ag. 4a is the mean lattice distortion of the solid 
solution. 


Figure 10 shows points for 4U in the series Ag—Zn, Ag—Ga, and Ag-Ge, 
all at 6atm®% of solute. For these, there is lattice contraction in the 
systems Ag—Zn and Ag—Ga, and a very small lattice expansion for Ag—Ge. 
The points lie on a straight line parallel to that for the Ag—Cd, In, Sn, Sb 
series but with a positive value of JU for da=0. These values suggest 
that, if we take a given atomic percentage of solute, there is in each 
series a similar linear relation between JU and da, but that the lines 
are displaced on passing from one series to the other, so that 4U increases 
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with the difference between solvent and solute in the Periodic Table. 
‘This relation is not, however, shown in the sequence Ag—Mg and Ag—Al 
where the former system has a large electrochemical factor as revealed 
by a maximum in the melting point of the Ag Mg 6 phase. ; 

As explained on p. 1186 the 4U terms are determined mainly by the 
depression of the % versus T curves in the systems concerned, and the 
above relations between JU and Ja are thus accompanied by similar 
relations between 4a and AT’, the depression of the z/7' curve below the 
melting point of the solvent. 


Fig. 10 
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AU is plotted against Aa for solid solutions containing 6 atomic% of Zn, 
Ga, Ge, and of Cd, In, Sn, and Sb. 


The relations described in the preceding sections confirm the view that 
in a series such as Ag(Cd, In, Sn, Sb) and Cu(Zn, Ga, Ge, As) the liquidus 
and solidus relations are determined primarily by electron concentration, 
but that this simple dependence is modified both by lattice distortion 
and electrochemical factor. The relative effects in the systems Ag(Cd, 
In, Sn) and Ag(Zn, Ga, Ge) suggest that, for a given atomic percentage 
of solute, the 4U values are a linear function of the lattice distortion 
Aa, but that each sequence involves its own characteristic type of 
distortion. In the systems Ag(Cd, In, Sn, Sb) and Cu(Zn, Ga) where 
the electrochemical factors are small, the correspondence between the 


we re 
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factors for the lattice distortions and depressions of the solidus curves is 


remarkable, and the apparently whole number relations require explana- 
tion. In a later paper it is hoped to deal with the relations in systems 


‘such as Cu-In where the ionic radius of the solute is large. 


ACKNOWLEDGMENTS 
Thanks are due to Professor Sir Cyril Hinshelwood, F.R.S., and to 
Dr. F. M. Brewer for laboratory accommodation and other facilities 
which have greatly encouraged the present work. Dr. A. Hellawell 
must also be thanked for determining the Ag—Ge liquidus and solidus 


‘points used in fig. 10. 


APPEND TX A 
THE SystEM CopPpER—ARSENICT 
1. Introduction and Previous Work 
The liquidus curve of the system Cu—As was investigated by Bengough 


-and Hill (1910), according to whom the «-liquidus falls steeply to a eutectic 


at 683°c, and then rises to a maximum in the region 25-28 atm% As, 


‘corresponding to the formation of the phase usually known as Cu,As. 
‘The «-solid solubility curve below 686°c was established accurately 


by the x-ray work of Mertz and Mathewson (1937) and of Owen and 
Morris (1949). In the present work the liquidus has been determined 


accurately almost as far as the second eutectic at 46 atm, As and 600°c. 


The «-solidus has been established by microscopical methods, and by 
x-ray methods some information has been gained of the solid phases 
in the region 25-32 atm, arsenic. 


2. Experimental Methods 


(a) Metals used 


The British Non-Ferrous Metals Research Association kindly presented 
three master alloys containing approximately 18, 36, and 46 atm% As. 
These were made from pure copper and high purity arsenic, and contained 
a trace of phosphorus in order to ensure freedom from oxygen. In view 
of the close resemblance between the chemical characteristics of phosphorus 
and arsenic, it was felt that the presence of a trace of phosphorus was 
less serious than that of oxygen. The Association also kindly prepared 
a series of cast and hot rolled copper-rich alloys which were used for 
the preparation of alloys of intermediate composition, and also for the 
solidus determination. 


(b) Cooling curve methods 

For alloys containing up to 30 atm, arsenic cooling curves were taken 
by conventional methods, using an atmosphere of purified argon in a 
closed tube furnace designed by Hellawell and Hume-Rothery (1957). 


+ The experimental work described was carried out by one of us (J.B.) for 
Part II of the Honours Course in Chemistry at the University of Oxford. 


P.M. 4M 
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The cooling curve ingots weighed about 40 g, and all the standard precau-- 
tions} were taken to ensure the accuracy of the temperature measure- 
ments and the absence of temperature gradients in the melt. The- 
ingots were weighed after the experiments, and up to 21 atm°% arsenic 
the loss of weight was negligible. Beyond this poimt the experiment 
was stopped as soon as the liquidus point had been established, and the- 
ingot was rapidly cooled so that no appreciable loss of arsenic occurred 
after the arrest. Separate experiments were made in order to determine: 
the temperatures of the eutectic and peritectic horizontals. All ingots. 
were analysed after the experiment. 

The above methods were satisfactory for alloys containing up to. 
about 30 atm%, arsenic, but beyond this point the loss of arsenic by 
volatilization was too great. Further progress was made by melting the 
alloys in closed silica tubes into which the thermocouple sheath was. 
inserted by means of a fused joint. The whole was packed in sand in 
an iron tube, and heated in an electric furnace which was rocked during 
the taking of the cooling curve. 


(c) The determination of the «-solidus 

The cast and hot-rolled specimens prepared by the British Non- 
Ferrous Metals Research Association were homogenized in sealed tubes. 
containing sufficient argon to give a pressure of approximately one 
atmosphere at the highest temperature used. After homogenization, 
specimens were annealed at successively higher temperatures and, after 
quenching were examined microscopically for the presence of chilled 
liquid. 


(d) Identification of phases by X-ray methods 

For this work use was made of the methods described by Owen and 
Morris (1949). Suitable proportions of the powdered master alloys. 
were heated together in sealed silica tubes at constant temperature, 
after which the tube was quenched, and a Debye-Scherrer powder x-ray 
diffraction photograph was taken. Owing to the volatility of the arsenic: 
a product of uniform composition was readily obtained after six hours 
at 575°C, and 24 hours at 450°c, and the phases present were identified 
by their diffraction patterns. This work was not carried out in detail, 
but it established the approximate positions of the phase-boundaries. 


(e) Analyses of the alloys 

The chemical analyses were carried out by Messrs. Johnson Matthey 
and Co. Ltd., and thanks are due to Dr. J. GC. Chaston for his care and 
attention. In the cooling curve work, the whole ingot was dissolved 
except for the small portions lost in preparing the microsections. For 


eS Ss SS ee eee eS 
+See W. Hume-Rothery and P. W. Reynolds (1937 a), and M etallurgical 

sedi Diagrams by W. Hume-Rothery, J. W. Christian and W. B. 
earson. 
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the solidus determinations the actual specimens used for microscopical 
examination were analysed. In the x-ray work the powders of the 
critical alloys were analysed after the annealing treatment at 575°o, 
and these compositions were assumed for the powders annealed at 450°c. 


3. Experimental Results 


The results of the cooling curve experiments are given in the table, 
and are plotted in fig. 11. The temperature of the first eutectic hori- 
zontal was determined as lying between 684° and 685°c from three 
cooling curve experiments, and as between 685° and 686°c from arrests 
on heating curves. The true value may be taken as 685°-+1°C, corres- 
ponding to the reaction Liquid=«+ 8. The liquidus rises to a maximum 
at 26-27 atm, As and 827°c at which temperature the B-phase is in 


Liquidus 
Atm% As point Notes Atm% As 
(c) 


1072 
1064 
1053 
1017 
969 
928 
860 
794 
735 
687 
757 


DWP Pere re 


A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 


Notes. The first and second, and fourth and fifth columns give respectively 
the compositions of the ingots by analysis, and the temperatures of the first 
arrests on the cooling curves. 

The experiments marked A were carried out by conventional methods 
and those marked B were made with sealed tubes (p. 1192). 


equilibrium with the liquid. The liquidus then falls to a peritectic 
horizontal at 709°c corresponding with the reaction Liquid+f=y, 
where the y-phase has the composition 30-31 atm As. From this 
point the liquidus falls to the second eutectic horizontal at 46 atm’, As 
and 600°c where the liquid is in equilibrium with solid y and solid arsenic. 

The «-solidus curve is shown in fig. 11, and was confirmed by two 
points from unpublished work by Miss Sutherland. The «-solid solu- 
bility curve in fig. 11 is that of Owen and Morris (1949) and shows an 
inflection at 380°C which suggests that a transformation on decomposition 
of the B-phase occurs at this temperature. No thermal arrests were 
detected in this region, suggesting that the transformation occurs too 


4M2 
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slowly or with too small an evolution of heat to be detected on an ordinary 
cooling curve. In agreement with the earlier workers an arrest was 
found at about 300°c for an alloy containing about 27 atm% As, but it 
was not possible to examine this in detail. 


Fig. 11 


Atomic As 


The x-ray diffraction films from the annealed and quenched powders 
gave clear and sharp lines for the 575°c specimens, and these suggested 
that the f-phase has a range of composition extending over at least 
latm%. At 450°c the films for alloys in the region of the B-phase were 
satisfactory and showed no significant difference from those at 575°c 
but for higher arsenic contents the films were less distinct, and it te 
hoped to investigate this problem later. 
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SouRCES OF Data 


For the liquidus and solidus curves, the tables of Hume-Rothery, 
Mabbott and Channel-Evans (1934) have been used, together with data 
from the following later sources. 


Ag-Cd, Ag-In, Ag-Sn, Ag-Sb. Liquidus Curves. Hume-Rothery 
and Reynolds (1937 a). 


Ag-Sb. Solidus. Hume-Rothery and Reynolds (1937 b) and (Miss) 
A. Sutherland (unpublished). 


Cu-Ge. Liquidus points by J. Reynolds (1955) agree with those 
of H-R, M, and C-E. 


Ag-Zn. Andrews, Davies, Hume-Rothery and Oswin (1941). 
Ag-—Ga.* Hume-Rothery and Andrews (1942). 

Ag-Ge. A. Hellawell. (Unpublished.) 

Ag-Mg.* Andrews and Hume-Rothery (1943). 

Ag-Hg. Smithells (1955). 

Cu-As. Liquidus and Solidus. Appendix A of present paper. 


Cu-In. Liquidus and Solidus. J. Reynolds (1955). The Liquidus 
is in good agreement with that of Weibke and Eggers (1934), 
but the depression of the solidus is much greater. 


The lattice distortion values have been taken from the data of Hume- 
Rothery, Lewin and P. W. Reynolds (1936), Owen and Roberts (1939), 
Owen and Edmunds (1938), and Owen and Rowlands (1940), and also 
from the sources marked * above. 
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SUMMARY 


An attempt is made to describe mathematically the state of stress at the 
‘centre of an edge dislocation. The discontinuity in stress across the slip 
‘plane suggests the use of the more general type of dislocation previously 
investigated by the author. As shown in the earlier paper the more 
general dislocation may be built up from an array of simple edge dis- 
locations. Here the converse is applied; the simple edge dislocation is 
presented as the difference of two or more dislocations of the more general 
type. By this means it is possible to eliminate the central singularity 
hitherto obtained in the elastic dislocation, to calculate the stresses at all 
‘points, and to plot the law of force acting across the slip plane. The 
method is believed to be new, and is found capable of wider application. 


§ 1. InTRODUCTION 


‘THE stress field of an edge dislocation as originally calculated by the 
theory of elasticity is not valid at the centre of the dislocation for two 
reasons, firstly because certain of the strain components are too large for 
elasticity theory to be applied, and secondly because the stresses and 
displacements as given by the theory become infinite at the centre. The 
problem is to find an alternative description for the state of stress at the 
centre consistent with the known elastic solution throughout the rest of the 
elastic body. 

Originally the singularity at the centre was avoided by cutting out a 
small circular cylinder around the dislocation line and applying elasticity 
theory only to the region outside this. But Orowan pointed out that the 
actual physical strains need be large only near the slip plane, so that the 
half space to either side of it might reasonably be treated as an elastic 
continuum. The result of this was the well-known calculation by 
Peierls (1940) later extended by Nabarro (1947) and by Foreman e¢ al. 
(1951) in which two elastic regions were joined across the slip plane 
by certain non elastic forces. These were related to the surface displace- 
ments of the two regions by a Law of Force, originally taken as sinusoidal 
in form but later modified by the insertion of parameters. 

By this method valuable results were obtained concerning the variation 
of width of a dislocation with the Law of Force, and the applied stress 
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necessary to set a dislocation in motion. But there are several funda— 
mental objections to the theory which cannot entirely be ignored. 

Firstly the solution does not agree with that of elastic theory even at 
large distances, although the strains are very small there and the theory 
might be expected to apply. Instead, Peierls’ displacements approach. 
those which would be obtained from the elastic solution by removing a 
strip of width a/(l—c) (where o is Poissons ratio, a the interatomic 
distance normal to the slip plane) along the slip plane, and rejoining the. 
remaining regions at a distance a while ignoring the differences in normal 
stress on the two surfaces. Thus even at great distances from the centre. 
the two half planes do not rejoin to form an elastic continuum, since even 
there the equations of elasticity are not satisfied across the slip plane. 

Secondly it is not always valid to assume, as Peierls has done, that the 
Law of Force must reduce to the form 


ee Vu 
a 


for small strains. Here zy is the shear stress or shearing force per unit 
area, . the shear modulus and Vw is the difference in the tangential. 
components of displacement : 


Ws ese Ab eee 


For in the elastic region the usual stress-strain relations must hold, i.e. 


The term du/dy corresponds to u Vu/a in the Law of Force, but év/dax has. 
been ignored. This term is by no means negligible for an edge dislocation, 
but is of the same order as du/dy along the slip plane. It is readily 
verified from the elastic solution (cf. Read 1953) that in fact, for 2 much 
greater than y, the relation between #7 and w is 


= - 2 Vu 
i 8 (2 On aS 
—2yu) a 


instead of the ‘ Law of Force’ as assumed by Peierls and Nabarro. 

The Law of Force was originally defined (Frenkel 1926) as the relation 
between applied stress and relative displacement when two half planes are. 
sheared uniformly over one another. Then no contradiction arises since 
there is no variation with x. The Law should not, however, be applied 
without modification to a non-uniform stress field. 

To avoid these difficulties a different approach is made in the present 
paper. The body is treated as elastic throughout but with discontinuities. 
of stress or displacement across a limited band of the slip plane. By 
assuming certain simple distributions of the relative displacement it is. 
possible to represent the centre in terms of higher order dislocations. The 
method of combining these dislocations to represent a Burgers edge. 
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dislocation at large distances is described in § 2. The method is applied 
to first order dislocations in § 3 and to second order dislocations in § 4. 
The resulting Law of Force is calculated and discussed in § 5. 


-§2. Toe GENERAL MeETHOD 


; As shown previously by the author (Mann, 1949) it is possible to find 
simple functions describing the state of stress in an elastic body having 
given discontinuities of displacement or transverse stress across a certain 
plane. A state of plane strain is assumed, and axes 2, y are chosen along 
and normal to the slip plane respectively as in figure 1. 


Fig. 1 


Diagram of an edge dislocation of arbitrary order showing coordinate axes. 


Let r, 6 be polar coordinates, measured from the positive a axis, and let 
u, v be the z-, y-components of displacement. Then an edge dislocation of 
order n is defined as one for which the relative displacement w’ is propor- 
tional to the nth power of the distance from the origin : 


“= [elee— [2 or 


et eee. oe ee ae Ss, AS) 
where «, is a constant. 
The Volterra or Burgers edge dislocation is described as one of zero order, 
since the relative displacement is constant over half the slip plane : 


TREE el a ee 


The stress system corresponding to an edge dislocation of any order n 
was derived in the previous paper. It may be expressed briefly in terms 
of a stress function U,, from which the components of stress and dis- 
placement are readily derived. 
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The method to be used in this paper is based on the following argument. 

Since all the equations of elasticity are linear, a stress or displacement 
derived from the sum of two stress functions is equal to the sum of the 
corresponding stresses or displacements derived from the stress functions 
separately. Moreover if all the functions are continuous and differentiable, 
it is possible to differentiate a stress function with respect to a, and, 
treating the derivative as a new stress function, derive from it new stresses 
and displacements. Each such component is equal to the derivative of 
the corresponding component derived from the original stress function. 

The stress functions of the edge dislocations afford particular examples 
of this property. The stress function U, of a second order edge disloca- 
tion is given by : 


apa 
Uy=By| 2e?—y") log r—dry 302+ al mane 46) 


where B is a real arbitrary constant, and from it are derived stresses and 
displacements : 
LYo=—4B(x log r—x—2y6) 
2. y=4BI(2—0)xy—(3—20)ay log r-+-6((2—c)y2—(1—o) x2). 
The relative displacement at 6=--7 is given by : 
2uu’ = — 87. B(1—o)x? 
47 B(1—o)a? 
LX 

=——Wof" SAV. 1S DAC) Bee Ser 


, 


uu = 


The function U, may be differentiated with respect to x to give a new 
stress function U,, 


© 0% 


=4By(x log r—x—y6) sw yin oe el en OF) 
which gives stresses and displacements : 


x y 
Ly = —4B( log r—- 2) 


am. 
2 4= —4B[(3—20)y log r+2(1—a)x 

aus 

RGD 
, Ou’, 
ETS. 

alee 87.B(1—o)x 
fe 
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U, is therefore the stress function of a first order edge dislocation. 
It may be differentiated to give a third stress function Ua 


U 
Uo = at 
ERTS Pua sear Ue ey) ee ee D) 
en Ow 
ee 
4Bx cos 26 
a ee ewer . . . . . . . (8) 


Ou 
2uu,= 2u— 
Mug 2. ox 
= — 2B[sin 20+-4(1—o)6]. 
The relative displacement is wu’, where : 


t du’, 
0 Oa 
paeosmb (laa) 
ea 
=2a%>. 


U, is therefore the stress function of a zero order or Volterra edge dis- 
location. Moreover comparison of eqns. (4), (6) and (8) shows that, 
with increasing order of the dislocation, the stresses and displacements 
become progressively more regular at the origin and less so at infinity. 

For this reason a more regular centre for the zero order dislocation 
may be formed by combining higher order dislocations at suitable positions. 
These must be so chosen that each component of the combined system 
approximates at large distances to the corresponding component of the 
zero order dislocation, that is, to its own first derivative if first order 
dislocations are used, or to its own second derivative if second order 
dislocations are used. 

Let some component of a first order dislocation take the value f(x, y) 
at any point P and values /;, f, at adjacent points P,, P, respectively 
arranged as in fig. 2. Then for sufficiently small h, and differentiable /, 
the finite difference relation holds : 

Ue are lis 

Ox Garris Ohm . . ° . . . . . (9) 
and, since f increases rapidly with r, (9) is also true for fixed h and suffi- 
ciently large 7, with proportional error of the order h?/r?. 


=f (2, Y2) 
where 25, 7 are the coordinates of P referred to D, as origin. Similarly 
fo=f(#1; 1) 


giving the approximate relation 


Of (x, : 
fle ye) Fen yx) LEY ee Oe ai) 
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for large r. This states that if a first order dislocation is situated at Der 
and a similar one of opposite sign at D,, then at any sufficiently distant 
point P a component of the combined stress systems is approximately 
equal to 2h times the derivative of that component, that is, to the corres- 
ponding component of the zero order dislocation. Such a system 1s. 
investigated more fully in § 3. 

Fig. 2 


Pater «hse a ee Side! 


D, O D ss 


<-———— 2k 


Diagram illustrating the equivalence between the difference at one point of 
similar components of two separated dislocations and the difference 
in the values at two separate points of a component of one dislocation. 


Similarly a representation in terms of second order dislocations may be 
obtained from the finite difference approximation to the second derivative. 
Second order dislocations are placed at D, and D,, and a double negative 
one at O. On dividing by h? the result at large distances is the second 
derivative of the component considered, that is, the corresponding 
component of the zero order dislocation (see § 4). 

The method may clearly be extended to third or higher order derivatives. 
if desired. Second order dislocations are sufficient for the present purpose 
since their components of stress and displacement are all finite at the 
centre. 

§ 3. First OrpER DIsLocATIONS 


The zero order edge dislocation with relative displacement (Burger’s 
vector) b may be represented by first order dislocations (5) by suitable 
choice of the arbitrary constant. Let one such dislocation be situated 
at D, and an equal negative one (B replaced by —B) at D, as in fig. 3. 
On dividing by 2h the stress function of the combined system is given by 

2B 
Us ee [(w-+-h) log re—(a—h) log r;—2h+-y(0,—05)] . . (11) 
where 7, 6, are measured from D,, and ry, 6, from D, as in fig. 3: 
r= (a+h)?+y? 


tan 0.= 25a 


etc. 
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From eqn. (6) the relative displacement takes the following values 
along the slip plane : 
u =0 forz>h 
bear, 47B(1—o)(h—2) 


ide TOD Wa go a (12) 
87B(1—c) 
| Oa ae for 2<—h. 
Fig. 3 


Be: 


Diagram showing the origins of the various coordinates. 


So u’=6b on the negative axis if B is chosen suitably : 
OM oa 

4 87(1—a) ° 

The stress components of this dislocation are readily derived from the 


stress function : 


Se | | 


Qn(l—o)| re rs 
a eee ls | ath «—h 
= | (a ra (13) 


ey —pb ry 9 1 1 
ee Onto) | los (A)+"( e 3) | | 


The stress components £% and 77 are finite everywhere, but 7 becomes 
infinite at D, and D,. The stress <% is discontinuous across the slip 
plane between D, and D, as indicated in fig. 4. Both displacements 
remain finite in the finite plane. 

At points a large distance from O the stress field approaches that of a 
zero order dislocation by (10, This may be verified directly. From 
fig. 3 

r, sin 0,=r sin 0=7, sin 0, 
r,cos@,th=r cos 0=rgcos0,—h . . . . . (14) 
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and, for r large compared with h, 


2h sin @ 
ioe (2) as 2h cos 6 
11 r 


with an error of the order /?/r2._ On substituting from (14) in eqn. (11) 
it is readily verified that U approaches U, for r large : 
2By | 2h cos 0 2h sin =| 

x: = 


Ur a , +2h log r—2h+-y. 


— Ors 


Similarly it may be verified that each stress component approaches that 
of the zero order dislocation at large distances. 
Fig. 4 


Curves drawn to an arbitrary scale indicating certain components of the 
edge dislocation expressed as the difference of two first order dislocations. 


The zero order dislocation is often represented by the symbol _|_, 
the horizontal line representing the slip plane and the vertical line the 
extra row of atoms above it. Ifthesymbol _, |is introduced to represent 
a first order dislocation as in fig. 1, the result of this section may be 
summarized diagrammatically as : 


wl] he ec le ae 
The final dislocation involves slip by the distance } if the first order 
ones have relative displacements -+-br/2h and are distant 2h apart. 
The symmetrical relation 

ees Pee 


is also valid. 
§ 4. SECOND ORDER DISLOCATIONS 


It was shown in § 2 that the stress function of the first order dislocation 
may be obtained by differentiating the stress function of the second order 
dislocation with respect to x. Therefore, repeating the argument of § 3, 
a first order dislocation may be represented at large distances by the 
difference of two second order dislocations. 

This may be stated symbolically : 


PIT che (aes |\ 0, Gee eee CED 


Centre of Dislocation : I 1205. 


where the symbol _11| represents a second order dislocation of the form 
shown in fig. 1. 


On comparing (16) and (15) a relation is obtained for the zero order 
dislocation in terms of four second order dislocations : 


ui) +] +] +a | =_L. Pe tendo) 
Positive dislocations of the type (3) are placed at D, and D,, and two 
similar negative dislocations at points nearer the origin, say at ~=-Lh. 
The relative displacement of the final dislocation is u’=b if the second 
order dislocations are such that 


2(h2— ke?) ° 

The simplest case, which will be considered in detail, is that of k=0, 
the two negative dislocations being situated at the origin. It will be 
shown in § 5 that varying k gives different Laws of Force. 

The result (17) may also be written in terms of stress functions : 


U=U, «+h, y)—U («+k, y)—U,(x—k, y)+ U(x—A, y) 
(12 12)U ((e, 9) 
or, for k=0 
U=U «+h, y)—2U,(x, y)+U(x—Ah, y) 
OO ae | mm tans Apter ers hat a Gh S| (18) 
In egn. (3) for U, the constant 6 must now be chosen as 
—pb/87(1—oa)h? 


Dt a| ev? log ((2-+8)?+¥?)—2e?—y") og (oy?) 


+ ((e—h)?—y?) log ((a—h)?-+y?)— 6h 


ae ord v4 =e 
—ty( (x48) tan deity tes 2x tan ae h) tan 4) | 


= on ba —y") )logr,—2y Ya,0,—32 |. oe een (9) 


The suffix 7 and summation sign have been introduced for brevity : 
Shi = fle +h)— fle) +f@—h) 


for any function f of x. 
The stress components are derived from eqn. (19) : 


en ee dO 9 | 
ee 27 (1—o)h? s| 8y ig (= 2)+2) 228 i 
ve —pb r 

CE he @n(1— oyna log ( af ; 


= sh a| Below —2y ¥ 4]. 


= 
bo 
= 

_—=— 
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At large distances these approach the stress components of the zero 
order dislocation (cf. Read 1953) : 


— —pb . 
Re eas 2 20 
ae ont 6(2-+-cos 26), 


= ub ; 
y= leg sin 6 cos 20, 


=~ pb 
Din d ain 26. 
2 bers cpargemmpars oe 6 


The components of displacement are given by : 
as b os 2h One 2 
v= pone Seem (1-+20)y? log ("1"2) +(1+20)h?+(1—2c) > a7 logr; 
87(1—o)h? 2 7 : 
+4oy > v8 | a), eo 
and approach those of the zero order dislocation at large distances : 


u os [sin 20+.4(1—o)6] 


Sal laa) 
The distributions of £z’, zy and wu’ along the slip plane are sketched 
in fig. 5. 


[cos 26+ 2(1—2o) log r]. 


er 


‘Curves drawn to an arbitrary scale for the edge dislocation expressed in terms 
of second order dislocations. The shear stress is finite and continuous 
all along the slip plane. 


On comparing this with fig. 4 it is seen that there is now no discontinuity 
in zx’, nor singularity in #y. The smooth curve of w’ is given by: 
u =) ore h 
b(h—a)? 


b(h?—2ha—a? 
ee! for —h<az<0 


= for =f, 


ie Bb vs li lO a 
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The maximum shear stress is about 1:766/27(1—o)h and occurs at 


r=+3V2h. This may be compared with the result of Foreman et al. 


(1951) 
+) ee ee 
RE Dex Mae ts 


where w is the width of the dislocation and approaches 0-6/ for large h. 


§5. Tor Law or Forcr 


From the values of u and £7 near the slip plane a curve may be drawn 
showing the relation between the relative displacement Vu of adjacent 
atoms and the shear stress acting at their centres. Peierls and Nabarro 
called this relation the Law of Force although, as shown in § 1, it is 
a law which changes form on changing the external stress field. 


Fig. 6 


The Law of Force for various dislocations. 


Carve 1 A=6-7b w=5d | 93) and (24 
Te 108 te (23) and (24) 
Lit h=6-76  w=0°7b 25) and (26 
eh 16h ap 12) oe (25) and (26). 
For comparison with previous work the Law of Force is again plotted 
in fig. 6. The ordinate is the value of xy on the slip plane. For con- 
venience in varying h the shear stress may be expressed as a function of 
ates 
oy. aly Ca (23) 


where ; 
= 1[(t—1) log ((—1)?—2# log #?+(¢+ 1) log (t+1)?]. 
Similarly Vu/a, for small a/h, is given by 


Vu ue (3—20)b a 9 
2. eee )(— ee a ere oe 
a a eta h Co 


P.M. 
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C= 0 tel 
=16(1—t)? Ox tasl 


= pols 


These formulae have been used in plotting curves I and II of fig. 6, 


with a=b, c=}. For other values of h, a and o the curves are readily 


obtained from the table. 


t OF 0°05 0229) O38 0-4 | 0-5 0-6 | 0-7 0-8 | 0-9 1 2 3 oO 


AM 0 | 0-66] 1-04]1-31 | 1-51]1-65 | 1-73] 1-76 | 1-74] 1-64 | 1-39] 0-54] 0-34] I/t 


u’/b| 0-5 | 0-40 | 0-32 | 0-245 | 0-18 | 0-125 | 0-08 | 0-042 | 0-02 | 0-005} 0 0 0 0 


Curves III and IV have been added to show how the form of the Law 
of Force may be varied by altering the positions of the second order 
dislocations. The two dislocations at z=0 are moved to z=%3h and 
x=—th respectively, and the constants are adjusted to keep the stress 
field unaltered at large distances. Then the shear stress on the slip 
plane is given by 


2ub 


Ce ior oy Ty (25) 


where 


2 ») 
x log (¢-+3)? + (t+ 1) log (¢+1)?] 


and 
VaR ua (83—20)b 
—=— 2 
% = Creag tie eter rh 2 ei) 
with 
2b 
CO ges Cost) 4<t<l 
2b 
=> (i-#) -<t<} 
2b P 
== (-2%-#) —1<t<—} 
=) aol 


Further variations may readily be obtained by choosing other positions 
for the central dislocations. In the limit as the inner dislocations 
approach those at =--A the solution reduces to that of two first order 
dislocations as in § 3. 

Narrower dislocations are obtained by decreasing h, but the width 
as measured along y=-+3b does not decrease indefinitely. Even for 
the limiting case h=0, when the stress field reduces to that of a zero 
order dislocation (Read loc. cit.) the displacement wu on y=+4b gives the 
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width 1-86. The width 1-5b used by Peierls is obtainable only by taking 
a=0-8b or less in the zero order dislocation. 

The elastic solution given in this paper cannot be expected to give 
so narrow a dislocation. The width 1-5) implies that the first row of 
atoms has been compressed or extended by the amount 15 within a 
distance 1-5), an average strain of 1/6, to which elasticity theory can 
scarcely apply. 


§ 6. THe STRAIN ENERGY 


The strain energy W per unit length of the dislocation is estimated 
from the work done on the sides of a cut from x=h to x=—R along the 
slip plane, assuming an outer boundary of large radius R. It is readily 
shown to be independent of the width h except for the usual term in 
log (R/h) 


h 
WaaeA | ey.w'de+Wr 
ae 


where Wp, represents the contribution from the outer surface and is 
independent of / if R is sufficiently large. From (23) and (24), 
R a 
xy.u' da +4] xy.u'da 
—h ee 
-— tu ra 
= —___. |" T. —t)T hd 
mins, T hdt+| (2t ) at 
pb? 
47(1—o) 
2 
fe a (toe +13), 


oh h 
=+| iay.n! dae =3 | 
—R 
pb? 
R 
= —_—— log a + terms independent of R, h 
~ 4ar(1—o) 


§ 7. CONCLUSION 


Having calculated the stress fields from various assumed distributions 
of relative displacement, it is desirable to compare the result with ex- 
periment and so choose the solution closest to reality. Unfortunately 
this seems difficult to do. The properties most apparent in experiment 
depend on long range effects, and at large distances the fields of all these 
dislocations are identical. 

The Law of Force at least provides a comparison with previous theories 
and has been calculated for that purpose. The curves I and II of fig. 6 
are similar in form to those of Foreman ef al., although with a different 
initial slope for reasons discussed in §1. Curves III and IV are more 
nearly of the form proposed by Orowan (1954). The maximum shear 
stress on the slip plane is of the same order as that obtained by Foreman 
et al. at the first row of atom centres. It is about yub/sh, the central 
singularity having entirely disappeared. 

The chief result of this calculation therefore appears to be a change 
in the initial slope of the force curve. 


4N2 
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But the method itself is of wider application and opens up a whole 
new field of problems of plane elasticity. Of particular interest are those 
concerned with inclusions, such as deposits of solute atoms along certain 
planes. If the dislocations of this paper are replaced by those of the 
second type where the discontinuity is v’ instead of wu’, the same 
method gives the stresses around an inclusion of very general shape. 
The problem of a narrow slit opening only very gradually from each end 
is now readily soluble, whereas previously it was treated as the limit of 
an ellipse with consequent difficulties and complications. The stress 
field of such a dilated slit is treated in a subsequent paper. 


ACKNOWLEDGMENTS 


I wish to thank Dr. P. B. Hirsch for very helpful discussion and criti- 
cism, Dr. J. 8. Koehler for reading the manuscript, and the Atomic 
Energy Research Establishment, Harwell, for a grant. 


REFERENCES 


Foreman, A. J., Jaswon, M. A., and Woop, J. K., 1951, Proc. phys. Soc. Lond. 
A, 64, 156. 

FRENKEL, J., 1926, Z. Phys., 37, 572. 

Mann, E. H., 1949, Proc. roy. Soc. A, 199, 376. 

Naparro, F. R. N., 1947, Proc. phys. Soc. Lond., 59, 256; 1952, Advanc. Phys., 
1352 

Orowan, E. O., 1954, Dislocations in Metals A.I.M.M.E., 74. 

PEIERLS, R. E., 1940, Proc. phys. Soc. Lond., 52, 34. 

Reap, W. T., 1953, Dislocations in Crystals (New York: McGraw-Hill), p. 116. 


i ae 


Photodisintegration of the Alpha Particle} 


By B. H. BRanspEN and A. GC. Doveras 
Department of Natural Philosophy, University of Glasgow 


and H. H. RopEertson 
CERN, European Organization for Nuclear Research, Geneva 


[Received June 21, 1957] 


ABSTRACT 


The cross section for the two particle, *He(y, p)?H, photodisintegration 
of the alpha particle is calculated, taking into account the interaction 
between the proton and triton in the final state. It is shown that this 
is insufficient to remove the discrepancy between the experimental 
cross section, and cross sections calculated with wave functions that are 
consistent with the binding energy of the alpha particle, even when the 
tensor force is included in the two body interaction. 


§ 1. INTRODUCTION 


Tue study of the photodisintegration of the lightest nuclei may, in 
principle, yield valuable information about the nuclear wave functions 
and interactions. One of the most suitable reactions for investigation 
is the two particle, (y, p) or (y,n) photodisintegration of *He, as the 
cross section has been measured in detail at low energies (Fuller 1954, 
Perry and Bame 1955, Atkinson et al. 1957). The cross section (fig. 2) 
rises rapidly from the threshold (at 19-8 Mev), to a maximum value 
of ~1-8x 1027 cm? at 6 to 7 Mev above threshold. The theory of the 
reaction has been developed by Flowers and Mandl (1951) and by Gunn 
and Irving (1951), who were able to account for the qualitative features 
of the observed cross section, without assuming the existence of an excited 
state of 4He. Gunn and Irving described the ejected proton by a plane 
wave and used simple analytic wave functions for *He and *H of the form : 
Gaussian 


$=N exp {—w(Y747)} ON (La 
i>j 
Irving : 
p=N exp {—p( > 1557) ( > 747)? Oe 10) 
i>j i>j 


They found that if the parameters » (denoted subsequently by py for 
the triton and py, for helium), were determined variationally to give the 
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best 3H and 4He binding energies, using central forces, the calculated 
(y, p) cross section exhibited a maximum at an energy much higher than 
that observed. To obtain a maximum at the correct energy, the size 
of the 4He nucleus had to be increased considerably (4, decreased), 
but when this was done the value of the calculated cross section at 
maximum became several times too large. 

In the present paper, the interaction in the final state is taken into 
account and it is found that the position of the calculated maximum is not 
significantly altered, although the detailed shape of the cross section 
does depend on the exchange nature of the nuclear force. Wave func- 
tions consistent with the binding energy of ‘He, for two body interactions 
containing a tensor component are used to recalculate the (y, p) cross 
section. It is found that the consistency problem remains and that it is 
impossible to choose wave functions, of this type, that give both the 
correct position of the maximum cross section and also the correct 
binding energy of *He. 


§ 2. GENERAL THEORY 


The electric dipole transition matrix element for the (y, p) reaction is 


M= >| Pier ee) Oars) eee 
spin 
and the total cross section o(y, p) is 
me? 
oly, P)=2 r) Kl Po sy eee ee 


where z,, 2, are the z components of the coordinates of the protons, 
referred to the centre of mass of the system, Y,, Y, are initial and final 
nuclear wave functions, m is the nucleon mass and K, k the photon and 
ejected proton wave numbers, respectively. 

By the conservation of energy, if the binding energies of 4He and 3H 
are denoted by ¢«,, ¢, then 


_ 2h? 
he = sh eq—€ g. 2b ty) age ee 


Assuming for the moment central forces, the initial “He wave function 
may be written 


W =4,(1234)o9(1234) Lies Tet) 


and the final (proton-+-triton) state may be represented by a resonating 
group wave function 


1 ; 
Py= Ty (1—Pra)fo(234)F (1)oo(1234) ee 7) 


where Po oy, are the ‘He and °H spatial wave functions, (1) describes 
the motion of the proton relative to the triton, Pi, is the exchange 
operator and the singlet spin function oy is 


Oo(1234)= 2 {a(1)B(2)—B(1)a(2) } {(3)B(4) —B(3)a(4)}. 


SA ee re NG 
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If ¢,, dp are symmetrical functions of the coordinates of all particles, 
then 


1 
M=—, | $,(r, X, R)dy(X, R)F(r)(r,+$X,)de dX dR. . (7) 


where the coordinates 
r=—"y+3(rotrstry); X=—reth(rg+ry), 
Rerirts have been introduced. As ¢,, dy represent S states, the proton 
is ejected into a p state, so that 
P(r)=(kr)-*f,(r)P1(cos 0) 


and 
fil(r)~sin (kr—da+a log (2kr)+58), 
and 
_ 2me* 
c= — 32k . 


§ 3. THe INTERACTION IN THE FINAL STATE 


Bransden et al. (1956) have calculated wave functions for neutron- 
triton scattering by the resonating group method assuming central two 
body interactions of either the Serber or symmetrical exchange type. 
This work is being extended to proton-triton scattering and radial 
functions /;(7) have been found over a range of proton energies up to 
40 Mev (Bransden and Robertson 1957). The triplet even two body 
interaction is Gaussian 

V(x)=A exp (—A?z?) Det ako) 
A=—45 Mev, )2=0-267 x 1026 cm-? 
and the triton is represented by the Gaussian function (1a), with 
px —0-0718> 102° em-2. 

Comparison of the calculated n—H* and n—He?® scattering cross 
sections with experiment suggests that (8) is a reasonable equivalent 
central interaction for the four body problem when a Serber exchange 
factor is used. 

The initial He ground state wave function is also taken to be of 
Gaussian form (1 a), and the constant , is chosen so that the maximum 
of the (y, p) cross section occurs at ~7 Mev above the threshold for a 
plane wave final state (f,=7j,(kr)). This gives »,2—0-021 x 106 cm-?, 
a value which leads to a very large root mean square radius, 
R=2-5x<10- em for the helium nucleus. 

In fig. 1, comparison is made between (y, p) cross sections calculated 
with plane wave final state wave functions and those calculated with 
wave functions /,(7) determined from Serber and symmetrical two body 
interactions, (8). It is seen that the interaction between proton and 
triton in the final state alters the shape of the (y, p) cross section 
appreciably, but that the position of the maximum given by the f,(7) 
determined from the Serber potential (which is consistent with the 


1214 B. H. Bransden, A. C. Douglas and H. H. Robertson on the 


n-t scattering data), is quite close to that given by the plane wave. The 
use of the symmetrical interaction gives a maximum at an energy about 
2 Mev higher. 


Fig. 1 
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Energy of incident y-ray above threshold (Mev) 


Calculated *He(y, p)?H cross sections for Gaussian wave functions: (a) with 
a plane wave final state, (b) interacting final state—Serber potential, 
(c) interacting final state—symmetrical potential. 


If w, is altered to the value that gives the best *He binding energy, for 
the potential (8) (u,2—=0-0789 x 102° cm-?, «,——26-7 Mev), then the 
position of the (vy, p) maximum cross section is much higher (~30 Mev 
above threshold). A resonating group ground state wave function of 
form (6) for 4He has also been investigated. This wave function gives a 
binding energy as good as the Gaussian function («,——26 Mev) but 
has a much better asymptotic form. However, in this case the maximum 
cross section occurs at an energy which is still too high (~20 Mev above 
threshold), for both plane wave and interacting final states. 

It must be concluded that the interaction in the final state is not 
responsible for the discrepancy between wave functions that are con- 
sistent with the binding energy of 4He, and those that are consistent 
with the (y, p) cross section. In investigating the consistency problem, 
it should be sufficient to employ plane wave final states and this has 
been done in the following sections. 


§ 4. Tur Tensor Force 


It has been suggested (Gunn and Irving 1951) that the discrepancy 
might be removed, if *He wave functions consistent with a two body 
interaction containing the tensor force (which is known to be important 
for the binding energy problem), were used in calculating the cross section 
In the presence of a tensor force, D state wave functions are directly 
coupled to the 5 state wave functions in both the 3H and *He ground 
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states. The electric dipole matrix elements (3) vanish for an S = D 
transition and the D = D transition is unimportant as the percentage 
D state of °H and ‘He is less than 46° and as the S = S transition is 
not abnormally small. Thus only the S parts of the 3H and 4He wave 
functions are considered and are used in the expression (7), which still 
holds. 

Irving (1953) and Abraham et al. (1955) have investigated the binding 
energies of 9H, *He with two body potentials containing central and 
tensor components, with Yukawa radial shapes. The S and D wave 
parts of the 9H and ‘He were represented by the exponential form 

o=N exp —p( > 1;,7)1? DAE y fen een rma 
>i 
where is a variational parameter. 
With this wave function, and with a plane wave final state 


e Sole 
oly, P)=3( Fa) { + Sys ler ea) 


Ua pug’ Riis eee 


a) F(z) 2 10 
x a oe 13 225 {f( )} ( ) 
oy 4: 16 
where 
1 I 
2)= — ——— (5852?— 9752+ 422 
I (2) 94 vote ve ae ) 
I 
— ——__ , —— (5 9422+ 66z+8 
Te seh a pet as any 
p 392*— 52z+ 16) sin-? Toke 
ae! z 2 sl AE 
and 


Poe +h?/3 
Pood h?/3— Bpuq?/4 
Abraham e¢ al. considered that if the nuclear potential is of the form 
exp (—Az) _ps,,°2 (—7x) anak 
AX Hx 


ae 


V (x)= —A[1+ 29(0, . ¢.—1)] 


then the set of constants, consistent with the binding energy of 7H, 
and with the low energy nucleon-nucleon scattering data, that is most 
likely to be also consistent with the *=H-and *He binding energies is 


A=39-51 Mev, B=33-03 Mev 
K-t=21°1847< 10-* cm, Heo loslo x 10-™ ent. 
The variationally determined parameter (including the principal D 
state only) is for this case wy=1-134Xx 101% cm™}, giving a_ binding 
energy <,——16-6mev. (This increases when further D states are 


included in the wave function.) The (y,p) cross section calculated 
from (10) with p y=py=1-134 x 104 em-! is shown in fig. 2, where it is 
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compared with cross sections calculated with w y= aus < 1033 em and 
[y=byp=0-:5 X 1048 cm and the experimental cross section. The cross sec- 
tion for u,—1-134 x 1018 cm rises much too slowly. If the tensor range is 
decreased slightly, the calculated binding energy can be raised to —28 MeV, 
but in this case j., is larger and the disagreement with the experimental 
cross section is worse. For ,—0-5x10%cm, the calculated cross 
section is in harmony with the experimental cross section, but this value 
of p,, is far outside any value consistent with *He binding energy. 


Fig. 2 
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Energy of incident y-ray above threshold (Mev) 
Calculated ‘*He(y, p)?H cross sections for exponential wave functions : 
(a) w,=1-1341013 cm! consistent with binding energy of ‘4He, 
(6) w,=0-8 x 1018 em-}, (¢) w,=0-5 x 1018 em}, (d) the experimental cross 


section (this is subject to an error of +20°%) taken from Perry and 
Bame (1955). 


§ 5. CONCLUSION 
To summarize the results of Gunn and Irving (1951) and of the present 
work, fig. 3 shows the variation of the position of the maximum cross 
section with the size of the He nucleus for the exponential (9) and the 
Gaussian and Irving (1a,b) wave functions. (The approximation 
/y=Hy has been made, but it should be noted that o(y, p) is relatively 
insensitive to jy.) It is clear that for wave functions that are functions 


of (>/7;;?) only, the position of the maximum depends sensitively on R, 
i>j ? 


the root mean square radius of ‘He, but is not very dependent on the 
shape of the wave function. 

The electron scattering experiments of McAllister and Hofstadter 
(1956) give a mean radius of 1-61 10-13 em for 4He. Wave functions 
that agree with this value are not consistent with either the calculated 
binding energy (which requires a lower value of R for the type of wave 
function considered) or with the calculated (y,p) maximum (which 
requires a larger R). The former result is true even for the more elaborate 


ee Oe 
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six parameter wave function of Clark (1954) (also a function of (> 777) 
where for a calculated binding energy of —29-87 mev, the mean radius 
is about 3 of the observed value. It seems likely that these inconsistencies 
may be resolved when hard core interactions are introduced into the 
nuclear potential, which tend to push out the wave function, for a given 
binding energy, thus increasing R and at same time decreasing the 
calculated (y, p) maximum}. Kikuta e¢ al. (1956) have considered the 


Fig. 3 
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Energy of maximum cross section above threshold (mev) 
The variation of the position of the (y, p) maximum cross section with nuclear 


size: (a) Gaussian wave functions, eqn. (1 a), (b) Exponential wave 
function, eqn. (19), (c) Irving wave functions, eqn. (1 6). 


effects of hard core interactions on the binding energies of 3H and #He 
and they find that because of the increased value of #& for these nuclei, 
the well known discrepancy between the coulomb and binding energies 
of ?He can be removed. The construction of *He wave functions from 
two body interactions containing a repulsive core is being undertaken 
and the photodisintegration of *He will be investigated with the new 
wave functions. 
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+ Since this work was completed, a paper by Rustgi and Levinger (1957) 
has appeared in which the total and bremsstrahlung weighted cross sections 
op=fo dE, o,=[(c/E)dE, are calculated for the Irving (1953) wave function. 
Although these cross sections relate to the sum of all dipole reactions in *He, 
they find that the o, inferred from the experimental (y, p) and (y, n) measure- 
ments is already considerably larger than that calculated, and they draw 
conclusions from this result similar to our own. 
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ABSTRACT 


By consideration of the observed spectrum of *Li in intermediate 
coupling, values are obtained for the p-shell radial integrals, the spin— 
orbit force parameter and two of the three independent parameters of 
the exchange character of the central force. These results are used to 
calculate the spectrum and properties of 7Li in intermediate coupling. 
Fresh evidence is found for the assignment of the 7-46 mev level in 7Li 
to the 5/2~ level of the *4P term, indicating the probable presence of a 
further 5/2~ level (from the ??F term) in the 5 to 7 Mev region. 


§ 1. InTRODUCTION 


THE success of the intermediate coupling shell model in the interpretation 
of spectra, magnetic moments and a wide variety of reaction properties of 
light nuclei leads us to consider the problem of obtaining, from available 
experimental data, information about a two-body central force mixture 
which, together with the well-known spin-orbit term of the model, will 
describe the spectra of such nuclei. Previous work on the model has 
almost invariably started with some assumption about the force mixture 
to be used and in the majority of cases this has been that of Rosenfeld 
(1948) as far as the central forces are concerned. 

We shall show in this paper that it is possible to derive from the observed 
spectrum of ®Li a number of conditions on the exchange parameters of the 
force mixture which reduce them to a single one which cannot be predicted 
from the available data. Values of the p-shell radial integrals Z and K 
are also obtained for this nucleus, together with the strength parameter a 
of the spin-orbit force. These results are applied to the structure prob- 
lems of *Li. 

Recently Meshkov and Ufford (1956) attacked a problem related to 
this by applying the method of Bacher and Goudsmit to the calculation 
of the elements of the energy matrix of 7Li from those of ®Li. The 
6L.i elements were estimated by assuming a central force mixture and 
using a successive approximation method to find the values of the radial 
integrals and the strength of the spin-orbit force which gave the best 
fit to the (over-determined) set of conditions obtained from the experi- 
mental data. These values were taken over by the Bacher—Goudsmit 
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method into the calculation of the 7Li levels. The condition that these 
should correspond to experimentally observed levels is then used to give a 
revised estimate of the central force mixture. 

Apart from its uncertainties, this approach suffers, as its authors point 
out, from the assumption that all the parameters involved remain un- 
altered from nucleus to nucleus. In view of the fact that other studies 
suggest that the intermediate coupling situation changes systematically 
from nucleus to nucleus through the p-shell, it would seem desirable to 
devise a method which is potentially free from this restriction. 

In the present paper we calculate all but one of the unknowns directly 
from the spectrum of *Li. The remaining parameter has a comparatively 
small effect on the low-energy spectrum of many nuclei, and an arbitrary 
choice is made for it before applying the results to “Li. We are then 
able if need be to treat the radial integrals and spin-orbit strength as 
free parameters in the spirit of the intermediate coupling model. We find 
in fact that a satisfactory fit to the known spectrum and properties of 7Li 
is possible using virtually the same values of the spin-orbit strength and 
the direct radial integral L, when expressed in units of K, as were found 
for the nucleus ®Li; only the value of K (which determines the overall 
scale of the predicted spectrum) needs to be altered appreciably. 

The derivation of the unknowns from the °Li spectrum is given in § 2 
of this paper. § 3 goes on to deal with the application of these results 
to “Li. 


§ 2. THE Forcr MIxTURE 


In this section we calculate a single-parameter force mixture to describe 
‘Li exactly in its low-energy features, on the assumption that these can be 
described in terms of a pure p? configuration. It may be hoped that such 
a mixture may prove valuable in the study of other light nuclei. We shall 
ignore any possible effects of the presence of tensor forces, since, although 
Elliott (1956) has shown that some tensor force is necessary to explain 
the beta-decay of /C on the shell model, there is no evidence to show that 
the strength of such a force is in any way comparable with the central and 
spin-orbit terms. The extra parameters which would have to be intro- 
duced in order to include this force might well lead to serious errors in the 
calculation, even if we could find sufficient experimental data to determine 
them formally. 

The data we shall use on the spectrum of ®Li are given in fig. 1 which is 
extracted from the compilation of Ajzenberg and Lauritsen (1955). 
The diagram also shows (where these are not known experimentally) 
the spin and parity assignments we shall make on the basis of the shell 
model for the first six levels. In L-S coupling they are the levels, aoe 
IST) ss 4 2 Sy, Do: 

Confining ourselves to central and spin-orbit forces, we see that we 
must consider six independent parameters. These are three out of the 
four parameters W, M, B, H of the central mixture (their sum is normal- 
ized to unity), the spin-orbit strength a which occurs in the expression 
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a>j\;-s; for this force, and the radial integrals L and K. The first six 
levels appear to give us five equations (for the spacings), although the 
breadth of the level at 5:4 Mev will introduce some uncertainty which 
we shall consider below ; a further equation may be obtained by computing 
the energy zero for the theoretical calculation from known data. This 
zero represents the energy of a proton and neutron in the presence of the 
central core, but in the absence of forces not due to this core. 


LON 

5:31 (1,2) ) ey 

WWI YI 
(0,2) 


Spectrum of *Li (after Ajzenberg and Lauritsen 1955). 


We denote this zero by °Liy. We then have the energy equation : 
SLi, —*He+p+n-+C(ns*)+C(ps*) 


where the last two terms represent the energy of interaction (including 
the Coulomb energy which, though we assume it to be common to all 
terms of a given nucleus, must be taken into account when relating 
different nuclei) of the p-shell neutron and proton with the core of four 
s-particles. To evaluate these terms we consider the observed resonances 
that represent the ground states of °He and °Li. These have the form 
of split 22P terms with very broad upper levels. We are interested in 
the positions of the ground states of these nuclei in the absence of this 
splitting, i.e. in : 


5Lig(++n)=*He+p(+n)+C(ps*) 
5He,(+p)=*He(+p)+n+C(ns*) 
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where we have added a nucleon on each side of each equation for the sake 
of comparison with SLi. As it is impossible to obtain the spin-orbit 
splitting in these nuclei with any accuracy directly from the resonances 
concerned, we assume it to have the value given directly by the splitting 
between the J —2, 3, 7’=0 levels in *Li (see below), namely 3q—=2-331 Mev. 
This gives us the *Li, and *He, levels lying at }a=0-777 Mev above the 
observed ground states, i.e. (with the added nucleon) at 6-277 and 
5-437 Mev respectively above ®Li ground. Also we have on this scale : 
*He+p+n=3-703 Mev. 
These equations show ®Li, to lie at 8-011 Mev above ®Li ground. The 
table lists the states of *Li, giving their energies measured (downwards) 
from this zero. 
All energies in Mev 


Energy 
from zero 


8-011 
5-822 
4-441 
3-491 
2-701 
2-611 


0 
0 
1 
0 
1 
0 


Rm bob Oo We 


We now consider the energy matrices themselves. These have been 
given, for example, by Elliott (1953), who uses parameters slightly different 
from ours. The simplest are the 7=0, J=2, 3 matrices which consist 
each of a single term : 


BD: L—K—4a; BD, : L—K-a. 
Equating these to the observed energies we obtain : 
a=1:554; DL=—K+4-268 


where here and subsequently all energies are expressed in Mev. 

With the remaining matrices we proceed by substituting the values 
above for a and L into the elements and subtracting from the diagonal 
elements the observed root given in the table. On writing down the 
conditions that the resulting determinants are zero we obtain four simul- 
taneous equations in AK and the three central force parameters, which can 
be solved. 

This procedure however is to some extent nullified by the great width 
of the level at 5-4 mev. Analysis of the equations involved shows} that, 
while the parameters a, L, K and the combinations : 


310—W+M—B—H 
8O—W-—-M+B—H 


; I am very grateful to Dr, J, P, Elliott for pointing this out to me, 
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(note 8O=W+ M+ B+ H=1 in this notation) are extremely insensitive 
to variations in the exact position assumed for this level, the combination : 


NQO—W—M—B+H 


shows extremely critical dependence on the position, and remains, for 
practical purposes, undetermined by the data. The values obtained for 
the remaining parameters are : 


a=1:554 Mev : 2430—0-46 
L=5:45 Mev : 8QO—0-15 
K==1-18 Mev 1s(4O=1),. 


These values are not consistent with those of the Rosenfeld (1948) 
mixture which is based on saturation conditions. It would seem prefer- 
able to use values based on experiment in light nuclei rather than those 
which are obtained by extrapolation from heavier nuclei where, indeed, 
saturation may be due to properties of the forces other than their exchange 
character. 


§ 3. APPLICATION TO “Li 


The difficulty of obtaining a value for the parameter 10 persists 
throughout the p-shell and into the (2s-1d) shell, as Elliott and Flowers 
have found in their work on the nuclei of masses 18 and 19 (Elliott and 
Flowers 1955) and on the odd parity states of 1®O and 1°N (Elliott and 
Flowers 1957). In the latter paper they do obtain a value 110 = 0-5, 
but this estimate is not directly applicable to our case as they use values 
of 340 and #80 which differ considerably from those we have found. The 
difficulty arises from the fact that the parameter 40 (which in the two- 
particle case refers to an antisymmetric space state) does not occur in 
most of the matrix elements of the highly symmetric wave functions that 
provide the main components of the low-energy states of these nuclei. 
Conversely it follows that many of the properties of these states will be 
largely independent of the value we choose for it. This is the case for 
example in 7Li. Meshkov and Ufford (1956) observed that variations 
of the assumed position of the 5-4 Mev level in ®Li did not make any 
appreciable difference to their *Li results, although, as we have seen, this 
assumption affects the value of 110 critically. The reason is that the 
11P(11) state of p?, which involves 140, is not a parent state for the 
22P(3), 22F(3) and 24P(21) terms in “Li which give rise to the low-energy 
spectrum in L-S coupling, and the resulting insensitivity will persist for 
some way into the intermediate coupling region. 

Bearing this in mind, we make the arbitrary choice ''O=0, which gives 
rise to the force mixture : 

W=0:40, M=0:33, B=0:17, H=0-10. 
It is interesting to note that this force from *Li is very similar to the 
approximate mixture obtained by Meshkov and Ufford from the spectrum 
of Li (W=M=0-4, B=0-2). We may thus expect it to give good 
results when applied to this nucleus. 


P.M. 7 see 
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In intermediate coupling the L-S states of 7Li listed above give rise 
to the level sequence 


F312) Wo Tie oy 20) 2 


in that order. The size of the 7Li energy matrices makes any ‘ fitting ’ 
procedure such as that described above prohibitively tedious, and a 
further deterrent is the fact that the experimental situation is by no means 
so clear-cut as in *Li. Accordingly we have adopted the procedure of 
fixing the value of the ratio L/K at the value (4-62) that we found for °Li, 
and treating a and K as free parameters. On general grounds we do 
not expect the value of this ratio to vary greatly, though it is possible 
that evidence of variation may be found in further studies of the system- 
atics of the p-shell. 

On writing the parameter a in units of K, the latter quantity appears 
only as a factor determining the overall energy scale. The matrices were 
set up in this form and latent roots and vectors extracted for ten values 
of the ratio a/K. This work was largely carried out on EDSAc, the Cam- 
bridge Mathematical Laboratory’s digital computer. 

Meshkov and Ufford have discussed the experimental situation in ‘Li 
in the light of shell-model theory and have reached the same conclusion 
as ourselves, namely that the 7-46 Mev level observed is the **P;/, level, 
and that the existence of a further 5/2 level between 4:61 and 7-46 Mev 
is indicated, possibly masked by the broad positive parity level at 6-6 Mev. 
With this assumption the spectrum was fitted to the theoretical curves, 
and the best overall fit was found to occur with a/K approximately equal 
to 1-4. This is so close to the value (1-32) found for ®Li that the two may 
be taken to be equal to within the accuracy of the calculation, though 
the increase may be a genuine reflection of that found by Inglis (1953) 
and Kurath (1956) to occur throughout the p-shell. The observed and 
predicted spectra are shown in fig. 2 where it is seen that a very satis- 
factory account of the spectrum may be given in this manner. The value 
of K necessary for this fit is 1-56 Mev, rather larger than that found for *Li, 

Some of the reaction and electromagnetic properties of these nuclei 
have been considered by Lane using Rosenfeld forces with L=6K. We 
have carried through similar calculations and many more with our 
present wave functions and find complete consistency with the parameter 
values given above except in calculations involving the usual form of the 
MI operator, i.e. magnetic moments and Ml transition widths. We hope 
to discuss these anomalies in a subsequent paper. 

One reaction property of this nucleus is however of considerable interest 
in that it supports the identification of the 5/2~ level at 7-46 Mev with 
that from the “P,,, term. Calculation of the reduced width of this 
level for nucleon emission to ®Li ground in intermediate coupling agrees 
with the experimental value for the 7-46 Mev level, while that of the level 
from “F,,, differs from this by a factor of three. The evidence for this 
identification thus appears to be strong on a number of grounds, 
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We thus see that a consistent and satisfactory account of these nuclei 
may be given in terms of the force parameters we have deduced. On 
going from *Li to "Li the effective strengths of the forces appear to 


Fig. 2 


SH IThe 52 


NANCE) 


AN WY 


PREDICTED OsseRvED 


Spectrum of 7Li. In the observed spectrum only states known with 
certainty are included. 


increase in much the same ratio, due probably to a modification in the 
radial shape of the wave functions, an effect we should expect to be in the 


direction observed. 
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ABSTRACT 


A procedure for dealing with localized defects in metals is proposed, 
based on three main assumptions. First, the validity of the Thomas— 
Fermi approximation is assumed, secondly departures of the perturbing 
potentials from spherical symmetry are neglected, and thirdly, no account 
is taken of any relaxation of the lattice which may occur on the introduc- 
tion of localized defects. However, the framework outlined is sufficiently 
general to enable the second assumption to be avoided, should this turn 
out to be necessary. 

It is expected that the method will be useful for foreign atoms entering 
the lattice substitutionally, provided the atomic size is such that no 
serious distortion of the lattice is produced, and also for vacancies, 
although in this latter case it will no doubt eventually prove necessary 
to consider relaxation effects. 

The method proposed is illustrated by studying numerically the case 
of a vacancy in copper, and an approximate self-consistent field is 
established for the perturbed lattice, based on the method of Slater and 
Krutter for the unperturbed crystal. A partial wave analysis using 
the free-electron approximation has then been carried out. The resistivity 
comes out considerably larger than that obtained by previous workers 
and possible reasons for this are discussed. 


§ 1. INTRODUCTION 


WHILE many problems remain to be solved in connection with the 
perfect periodic lattice, a good deal of attention is now being devoted 
to the perturbations introduced in originally perfect crystals by foreign 
atoms, interstitials, vacancies and dislocations. Particularly with the 
current importance of the problem of assessing the effects of irradiation 
on metals, the theoretical study of vacancies and interstitials especially 
seems a matter of some urgency. The present investigation represent 
an attempt to devize a somewhat general method of approximating 
to the self-consistent field and the electron distribution in metals in 
which localized defects have been introduced. Thus, dislocations are 
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immediately excluded from the present discussion. Furthermore, and 
of more immediate concern, it will emerge that the assumptions which 
at this stage it turned out to be highly desirable, if not entirely essential, 
to make, probably also exclude interstitials, or rather preclude the 
possibility of obtaining a realistic description of the perturbed crystal 
in this case. However, a modification of the present scheme appears 
possible for this case, and we hope to deal with this later. Thus, we 
think that the present approach will provide a first approximation for 
foreign atoms entering substitutionally in the lattice, and for vacancies. 

It is appropriate first of all to review the most important attempts 
which have so far been made to deal with the problems of concern to 
us here. Most of these have been based on the treatment introduced 
by Mott (1936), in which only the conduction electrons are considered, 
the Sommerfeld free-electron model being used to describe the unperturbed 
metal. The effect of foreign atoms in a monovalent metal say, was then 
described by introducing a point singularity with charge Ze at the 
position of the foreign atom, the impurity atom having (Z+ 1) electrons 
outside a closed shell. The redistribution of the conduction electrons 
and the self-consistent perturbing field were then determined by the 
Thomas—Fermi method, and this approach has been the basis of a good 
deal of more recent work. A detailed review has been given by Friedel 
(1954), and for the most important developments since this was written 
the article by March (1957) may be consulted. 

While such an approach has a good deal of value, particularly from 
the point of view of obtaining physical insight into the problems with a 
relatively simple model, it is apparent that there are cases with which, 
by its very nature, it cannot deal. Thus, such an approach, without 
modification, cannot explain why small concentrations of gold atoms in a 
silver lattice lead to an increase in the resistivity. Furthermore, when 
dealing with vacancies in other than metals formed of light atoms, any” 
incorporation in the valence scheme inevitably involves approximations - 
which must be rather drastic, particularly concerning the form of the 
perturbing potential near the vacancy. Such a situation, in fact, also 
obtains when the valence treatment is used for substitutional impurities 
which are not immediately adjacent to the solvent atom in the periodic 
table. 

Previous workers have made attempts of various kinds to overcome 
these limitations, and brief reference to their work will now be made. 
Thus Huang (1948) in an important early paper, attempted to calculate 
the excess resistivity due to a monovalent atom in a monovalent metal. 
In this work also, use was made of the Thomas—Fermi approximation, 
but in a rather less fundamental way than in the present paper. In 
connection with the objections to dealing with a vacancy in the valence 
scheme above, attempts have been made by Jongenburger (1953), and 
more recently by Blatt (1955 a) to improve the approximation to the 
perturbing potential by using the negative of the Hartree potential for a 
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positive ion for copper, and adjusting this in order to satisfy the sum 
rule of Friedel -(1952). These treatments appeared to lead to quite 
reasonable results for the excess resistivity due to vacancies (about 
1-3-1-5 micro-ohm cm peratomic °%), but the fields obtained in this manner 
can hardly be justified in a fundamental way. Further, although these 
authors concluded that the exact form of the perturbing field did not 
appear to matter very much in the calculation of the excess resistivity 
(Friedel 1954), Blatt (1955 b) has shown it to be of some importance when 
changes in thermeolectric power due to vacancies and interstitials are 
calculated. 

Finally, we should mention some unpublished work of Howell, which 
was kindly made available to us, and undoubtedly has considerably 
influenced the course of the present study. Howell attempted to deal 
with the vacancy problem from the point of view of calculating the energy 
needed to create a vacancy in copper. He concluded that the method 
he proposed, based on the Thomas—Fermi approximation, was not 
satisfactory for the calculation of this quantity. However, we think 
that the present work circumvents the main difficulties associated with 
Howell’s approach, although we focus attention on the perturbing potential 
and do not attempt energy calculations. 


§ 2. THe PEeRFecT LATTICE 


For the perfect metallic lattice we have to solve the fundamental 
Thomas—Fermi (TF) equation 


V2(V — V.)=n(V — V,)?/2 SRE Crs eee os fey) 
where 
32772 
PS 3h (2m)*/2, 


subject to appropriate boundary conditions. Hence V is the electro- 
static potential and —eV, is the Fermi energy. For a metal in which we 
construct the usual atomic polyhedra surrounding the atoms, the 
problem of finding an exact solution of (1) is a formidable task and has 
not as yet been solved for any lattice. However, as Slater and Krutter 
(1935) emphasized in their original paper on the application of the TF 
method to metals, in many cases the atomic polyhedron is of high 
symmetry and may, to a good approximation, be replaced by a sphere of 
equal volume. It is then consistent to suppose that within each atomic 
sphere the potential V is spherically symmetrical, and the condition that 
V is given by (1) is periodic with the period of the lattice is replaced by 
the simpler requirement that 
dV 
(22) or 


where & is the radius of the atomic sphere. This also ensures that the 
sphere is electrically neutral. 


Ww. 
‘ 
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At the centre of the sphere, where there is a nucleus of charge Ze, 
we have 


Ze 
V—V)>—as r+ 0, Pema me 1S ark eet 3 | 


In this case, eqn. (1) reduces to the ordinary dimensionless TF equation 


dp 3/2 
where 
Ze 
V— yo= — d, : . . . : . . - (5) 
Saal eae ee SMe eer eNews tr kid) 
with 


4 2 \2/3 h2 
=(s3) 2me2Z Us ° 


The boundary conditions are now 


$(0)=1 Weg ko eas, Sa ee Ek 
P(X) dd 
Jn Nobeeeete eee 


where 2, is the radius of the atomic sphere in dimensionless units. This 
type of solution has been made the basis of the present treatment of 
defects in metals; should it eventually prove to be necessary to go 
beyond the sphere approximation for the perfect lattice, the method 
described by March (1957) may be used to calculate the effect of angular 
terms in the potential. 


§3. THE PERTURBED LATTICE 
In this case, suppose a single defect is introduced at a lattice site 
which we take as origin of coordinates, and let the nuclear charge change 
from Ze to Z’e at this point. The new TF equation to be solved is 


V4(V+V,—V.)=4H(V+V,—S),)3” eran elo} 
where V,, is the change in electrostatic potential caused by the introduc- 
tion of the defect. Then we must solve (9) subject to the new boundary 
conditions 


/ 


Ze 
V+V,>—asr— 0, Ee ee aed) 


V ,, > 0 at infinity. ee ee LY 


The neighbouring atoms are assumed unaffected by the introduction 
of the defect ; that is, any relaxation of the lattice is neglected in this 
approximation. We can therefore focus attention on an atomic sphere 
of radius R centred on the origin, and within this sphere we should 
expect that both V,, and V are spherically symmetrical to a good approxi- 
mation. Then (9) again reduces to the dimensionless 'TF equation, 
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although the boundary conditions are now different. However, the 
boundary condition at infinity cannot be applied without examining 
the solution outside the central sphere and we shall now proceed to do 
this. 

Since we are then rather far from the defect we may assume 


Veo dee Log) ee ey 


and this must eventually be checked in practical applications. Then 
we may write (9) in the form 


V2V+V,—V.)=u(V—V_)?? E 5 (=) +0 ({ vr} ) | , (13) 


or using (1) 


Viv =tulV—V,)eV ee ee 


To solve (14), we now imagine V, and (V—V,)!? to be expanded in 
spherical harmonics about the origin. Of course, (V—V ))'/? is not a 
spherically symmetric function when it is expanded about any other 
origin than the centre of the atomic sphere in which the point under 
consideration lies. However, we assume that the angular terms in V, 
are unimportant, and then if the radial term is indeed a sufficient 
approximation, we obtain this by solving the equation 


La pee eee 
saa Vs)= 5h Va) ee 


where (V—V,)!/? is the average of (V—V,)1/? over the angular variables, 
and this may be obtained numerically as discussed in § 4. 

It will be seen that once the solution of (15) which tends to zero at 
infinity is known, the remaining constant is a simple multiplier, and this 
greatly simplifies the problem. It means that only one solution of (15) 
is needed, and the arbitrary constant is then obtained by matching at R 
with the solution for V,, in the central sphere. The way this works out in 
practice for a particular case is discussed in §5. 


§4. DETERMINATION OF (V—J,)}/2 


(V—V,)?? is given as a function of r by 
2 pm 
| (V—V,)¥® sin 6 dO de 
0 0 
2 pr 
| sin 0 dO dé 
0 0 
where r, @ and ¢ are spherical polar coordinates. This may be evaluated 


numerically for values of r >R by considering a network of points on the 
surface of a sphere of radius r with the origin as centre. Each of these 


(16) 


ee eee eee 
+t The following treatment is quite different from that of Howell and, we 
believe, avoids the difficulties associated with his approach. 


Fields round Localized Defects in Metals 1231 


points may be associated with the appropriate atomic sphere, and the 
corresponding value of (_V—V,)/? may be obtained by the method of 
Slater and Kruttert. The integrations in (16) may thus be carried 
out numerically for spheres of increasing radii, although it is necessary 
to employ a series solution for (_V—V,)!/* in the regions very near to 
atomic nuclei. For r<R, (V—V,)"? is the same as (V—V,)!/2 in the 
present approximation. 


Fig. 1 


rin d.u, 


(V—V,)* as a function of r (atomic units used), 


The results thus obtained for copper are displayed in fig. 1. The 
radius R of the atomic sphere corresponds to the point of zero derivative 
of (V—V,)*, the discontinuity in the derivative occurs at the nearest 
neighbour distance and the curve is drawn out to the next nearest neigh- 
bour distance. In view of the sphere approximation employed to obtain 
(V—V,)1/?, it was found necessary to make a slight adjustment of the 
curve between FR and the nearest neighbour distancet in order to match 
satisfactorily with the results in the central sphere. The correct lattice 
structure (face-centred cubic) enters the calculation of (V—V,)!/?, 
although it probably has a minor effect on the form of the final curve. 


§ 5. CASE OF A VACANCY IN COPPER 
5.1. Calculation of Perturbing Potential 


We shall now illustrate the general method outlined in §2-4 by 
describing a calculation we have made for vacancies in a copper lattice. 
As remarked earlier, relaxation of the lattice, which undoubtedly occurs 


+ A solution of eqn. (4) corresponding to a value of &=2-66 atomic units was 
kindly made available to us by Dr. R. A. Ballinger. 


¢ At no point was V—V,)!/? changed by more than 3 or 4%. 
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when an atom is removed, is neglected in the present work; some 
discussion of the errors which are thereby involved will be given later. 

In the central sphere, where V+ V, is assumed spherically symmetrical, 
(9) may be transformed to the usual dimensionless form (4) by putting 


(V+V,—Vo= (a), < Hot ey tele 


where x is defined by (6). To give the required behaviour of V, near 
the origin, the boundary condition for y is 
y(0)=0. of Ob, beeen, ee on 
Near the origin the following series solution may be used7{ 
y= Az[1-+b507-+0,24-+ 070% ... | re ke) 


where the coefficients are given in table 1, and A? is the slope of y at the 
origin. 


Table 1 

4. pyle sos 

3G 11” 14256000 

ZA vee 2719 A® 
580) 13°” 53374464000 
ae AB ee 2062147 

7 1440 15°” 11208637440000 

4 
31104 


The convergence of this series is good; for example for 4 ~0-06 it 
can be used out to x=8. 

Next eqn. (14) is solved outside the atomic sphere. The boundary 
condition that V,, shall be zero at infinity was approximated by putting 
it zero at the next nearest neighbour distance. The averaging described 
in §4 was thus needed out to this point. The function (V—V,)!” has 
a discontinuous derivative at the values of r equivalent to nearest and 
next nearest neighbour distances, due to the presence of the atomic nuclei. 
This complicates the solution of (14) somewhat, but the difficulty was 
overcome by first solving it in the region between the nearest and next 
nearest neighbours by the relaxation method (see, for example, Fox 1949) 
with the boundary condition that V, was zero at the next nearest 


+ This solution was obtained in connection with the present blem b 
Howell, and independently, as far as 6,3, by one of us (N-HM,) during ine 
course of earlier unpublished work on the molecules P, and AS;. 
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neighbour point. A Taylor series was then developed around the nearest 
neighbour point and the integration extended inwards to the surface of 
the central atomic sphere by a method due to Fox and Goodwin (1949). 
The solution was thus determined for the outer region apart from a 
multiplying constant. 

In practice the solution of (9) inside the atomic sphere was tabulated 
for a chosen value of the slope of y at the origin, using the series given in 
(19). This solution can then be transformed very simply to obtain the 
correct matching of solutions across the surface of the central atomic 
sphere. This is possible since, if y(a) is a solution of the dimensionless 
TF equation, then it is easily shown that ¥,(x,) is also a solution provided 


yy=ay | 
2 —de : ee ei SU) 
ad?=1 j 


By suitable choice of a, matching of the two solutions was achieved. The 
perturbing potential due to a vacancy is then determined uniquely and 
is shown as a function of in fig. 2. 


Fig. 2 
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Magnitude of perturbing potential as a function of r (atomic units used). 


5.2. Electron Distribution round Vacancy 
We also thought it of interest to examine the electron distribution 


in the atomic sphere surrounding the vacancy. In atomic units, the 
electron density in the perturbed lattice is given by 
93/2 
n(r)— 3,alV + eeeivnl tape oh oe 4 (20) 


Using the results obtained here for V,, and the Slater and Krutter values 
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of V—V,, we find that there are very nearly two electrons inside the 
central atomic sphere. Thus, on removing an atom, there is a significant 
movement of electrons into the ‘vacancy sphere’ according to the 
present model (see, however, the remarks in §5.3 below). 


5.3. Comments on Approximations Involved 


In the vacancy case, the form of the perturbing potential is such that 
the application of the unmodified Thomas—Fermi approximation seems 
open to some question. The introduction of a correction along the lines 
proposed by von Weizsicker (1935) therefore merits consideration, but 
the problem is then much more formidable than that considered here. 
In view of this, and bearing in mind that there is a degree of uncertainty 
with regard to (a) the magnitude of the angular terms in the perturbing 
potential and (b) the effect of relaxation of the lattice, we have not 
attempted to make quantitative calculations of such a correction. 
However, in the light of previous work using the von Weizsacker method, 
it seems that while the results for the electron density (which depends on 
the second derivative of the potential) discussed in § 5.2 above may be 
rather significantly modified, the perturbing potential, which is our main 
concern here, is unlikely to be markedly changed. If energy calculations 
were to be made eventually on the basis of the present treatment, much 
more careful study of the von Weizsacker correction would seem to be 
called for. 

A less fundamental approximation, the linearization employed outside 
the central atomic sphere, may be checked for the vacancy case using 
our numerical values for V,, and at the boundary of the central sphere 
we find V,/(V—V,)=0-19. The first neglected term is (V,/{V—V,})? 
in eqn. (13) and it is soon verified that this approximation is not serious. 


§ 6. CONNECTION witH PREVIOUS WoRK ON VACANCIES IN COPPER 


Particularly with a view to establishing some relationship between the 
perturbing potential we have obtained and the forms used earlier by 
Jongenburger and Blatt, we have made calculations by the partial wave 
method of the scattering cross section due to V,. Following these 
workers, we have adopted the free-electron model for the conduction 
electrons in calculating this cross section. However, it should certainly 
be borne in mind that for vacancies in copper the perturbing potential is 
so strong and has such a marked effect on the electronic wave functions 
near the vacant site that it may turn out to be a drastic approximation 
to compute the excess resistivity by treating the conduction electrons 
as free. Thus the phase shifts calculated by this means may differ 
significantly from those which would result when the same perturbing 
potential V,, was used in conjunction with more realistic wave functions 
for the conduction electrons. The present calculations should therefore 
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be regarded as semi-quantitative only, but some interesting points emerge 
from a comparison with previous work and we shall therefore briefly 
report our findings. 

The first four phase shifts corresponding to /=0,1,2 and 3 were 
determined and these are given in table 2, together with the value 
obtained for the excess resistivity for 1 atomic °/ of vacancies. This 
is considerably larger than the values obtained by the earlier workers 
(in the range 1-3-1-5 micro-ohm cm per atomic °%), and we shall now 
attempt to analyze the possible reasons for this. We note first that 
previous workers have made no attempt to obtain the perturbing fields 
from first principles, but invoked the Friedel sum rule as an aid in 
constructing the potential. This may be stated (Friedel 1952) in the 
form 


) CO 


Ware ©) (2a ly aa te eeihs vet et (22) 


af 10 
where 7, is the phase shift for the partial wave order / for electrons at the 
Fermi surface, N being the ‘screening charge’. For the case of a vacancy 
in a monovalent metal, N has been taken as —1 by previous workers. 
However, in the present calculation the perturbing potential is completely 
determined by our approximate self-consistent field calculation, and 
using the phase shifts given in table 2 we find 


N=— 2-25. 


Whereas previously the conclusion had been reached that the form 


Table 2. Phase Shifts and Excess Resistivity 
U 0 1 2 3 
m Pai S 20 0-44 eee of) P22 20-04 


A,=3-09 micro-ohm cm/atomic %. 


of the perturbing potential was not very important for the calculation 
of the excess resistivity, provided the phase shifts satisfy the sum rule 
(see, however, Fujiwara 1955, 1956), it is not surprising that we find a 
considerably greater increase in resistivity with a potential which leads 
“to N=—2-25. We feel therefore that until the status of the sum rule 
is clarified, some caution should be exercised in using it as a means of 
determining perturbing fields. 

Finally, we must attempt to assess briefly the possible effect of allowing 
for relaxation of the lattice. Here there appears to be some measure of 
disagreement between various writers. For example, Roth (1957) 
suggests that the work of Overhauser and Gorman (1956) overestimated 
the effect of relaxation by an order of magnitude. From some rough 
estimates we have made, we do not feel that the resistivity as calculated 
here will change by more than 0-5 micro-ohm cm per atomic ° due to 
relaxation effects. It seems highly probable therefore that the use of the 
free-electron approximation for the calculation of the scattering cross 
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section is much more serious than the neglect of relaxation, and is causing 
the excess resistivity to come out too high (see especially the discussions of 
Blatt (1955 a) and Overhauser and Gorman (1956)). A calculation of the 
excess resistivity using the present perturbing potential together with 
a more realistic description of the conduction electrons would be of very 
considerable interest. 
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ABSTRACT 


Rods of annealed polycrystalline copper (99-991°%) were deformed in 
tension at 77, 195 and 293°K. The initial mode of hardening was linear 
up to a temperature-depended tensile stress o,,; the mechanism was the 
same as the one leading to linear hardening (stage IT) in single crystals. 
The assumption that at stresses above a, ‘ yielding’ took place by the 
spread of slip from some of the small slip zones characteristic of the 
linear stage, and that this resulted in interaction hardening, is shown to 
account satisfactorily for the parabolic shape of the stress-strain curves 
above ov. 


§ 1. INTRODUCTION 


In a recent paper (Feltham and Meakin 1957 a), which will be referred 
to as I, the present authors showed that a linear stage of hardeningt 
preceeds parabolic hardening in annealed polycrystalline copper deformed 
in tension at room temperature. This hardening was interpreted in 
terms of the model proposed by Friedel (1955), while hardening in stage 
III was thought to be of the interaction type, and to arise from dislocation 
pile-ups as formed above a characteristic temperature-dependent stress 
oy at which Cottrell-Lomer barriers ceased to be effective in confining 
slip to the small rectangular slip-zones active in stage II. 

The primary object of the present work was to study work hardening 
in copper under similar experimental conditions but at lower tempera- 
tures, and to examine in how far hardening in stage III could be 
accounted for in terms of a simple model based on the interaction 
mechanism. The choice of low temperatures was made in an endeavour to 
minimize interference from relaxation and creep phenomena ; for similar 
reasons the tensile strain « was in no case allowed to exceed about 10%. 


§ 2. MATERIAL AND METHOD 


Polycrystalline rods drawn from oxygen-free high-conductivity copper 
(99-991%), which were used in the experiments, had a diameter of 2 mm 


+Communicated by the Authors. 

+This linear hardening is now generally referred to as stage II, to distinguish 
it from linear hardening due to easy-glide : stage I. The parabolic hardening 
discussed in this paper will be referred to as stage IT. 
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and were 10 cm long. They came from the same stock as the rods used 
in I. After annealing in vacuo for 2 hours at 600°C and furnace cooling, 
the rods were extended in liquid nitrogen, powdered carbon dioxide, or 
in air, by means of a motor-driven Hounsfield Tensometer at a constant 
tensile strain rate of 10-3 sec-!. The spatial mean grain diameter D was 
2-84 10-3 em. A more detailed account of the experimental procedure 
has been given previously (1). 


§ 3. EXPERIMENTAL RESULTS 
The stress-strain curves obtained at 77, 195 and 293°K are shown in 
the figure. The strains up to the tensile yield stress cy were too small 
to be measured reliably, and the measured points have not been recorded 
in the figure. The curves were found to be reproducible to a high degree 
of accuracy. 


o KG/CM «IO 


2 ] 2 4 5 6 7 

Ele 

True stress-strain curves of annealed polycrystalline copper. Full circles - 
represent the parabola given by eqn. (1). 


A linear part, terminating at a tensile stress og indicated by arrows 
in the figure, is seen to extend up to tensile strains of about 19, incall 
cases. This part of the curve is then followed by a parabolic section. 
The origins (o9, €9) of the parabolae and the coefficients of work hardening 
X» were derived by analysing the curves by means of the equation 


(5 — o9)?= x 9(e—<9) 
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as in I. Again, the parabolae seemed to have their origins close to the 
origin of the coordinate system ; the parabola 


o2=3-62 x 107(e—5 x 10-4) (kg/em?)? ee at) 


for example is shown fitted to the lowest curve in the figure. The values 
of oy, og, og/oy and of the shear modulus G@ are given in table 1. @ is 
the unrelaxed modulus, the temperature-dependence of which was 
assumed to be equal to that of Young’s modulus at low temperatures as 
given by Lozinskii (1956). 


Table 1 
Ah o OCG og/oy Gx 10-5 
re kg/em? kg/cm? kg/cm? 
77 384 800 2-09 5-18 
195 351 670 1-91 5-00 
293 313 550 1-76 4-83 


§ 4. THEORY 


4.1. The Yield Stress oy 


On the basis of the theoretical treatments of the relations between 
the parameters characterizing yielding and work hardening in single 
crystals and polycrystals respectively (e.g. Sachs 1928, Taylor 1938, 
Kochendorfer 1947, 1950, Bishop and Hill 1951) the tensile yield stress 
oy should be given by 

oy=cT,+S, Se ene te aah A) 


where 7; is the critical resolved shear stress of the single crystal, and 
2<c<3. No adequate explanation of the significance of the stress S 
exists. Kochendérfer sought to relate it to elastic deformations due to 
hydrostatic stresses near grain boundaries. He shows that such an 
interpretation associates S with the Bauschinger effect, i.e. a polycrystal 
which had been deformed in tension just above oy should subsequently 
have a yield stress in compression given numerically by 


(oy )soRp SOT pm ee 


Such an interpretation of the significance is however untenable in the 
present case, for oy is about ten times greater than 7, as obtained by 
Blewitt et al. (1955) with single crystals of copper of purity comparable 
with that of our polycrystals, so that S is numerically about three times 
greater than cr;, leading to the inadmissible result of a negative numerical 
value for (cy)comp: AS it is reasonable to assume that stresses due, 
for example, to cooling from the annealing temperatures, are induced 
near grain boundaries in an elastically highly anisotropic metal such as 


P.M. is 
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copper, a small plastic deformation, possibly involving slip on several 
slip systems in each grain and the formation of some Cottrell-Lomer 
barriers, may have occurred. Such a hypothesis would not be open to 
the same objections as Kochendorfers ; a detailed theory of the structure 
of grains in annealed polycrystals has not however been developed as 
yet. 
4.2. Linear Hardening 

The small temperature-dependence of the coefficient of linear hardening 
¥(=do/de);, apparent from table 2, supports the earlier conclusion (I) 
that hardening in stage II arose through the locking of slip by Cottrell— 
Lomer barriers (Friedel 1955, Seeger et al. 1957). Friedel pointed out 
that the product o,7'/? should vary but little with temperature if og 
is the stress at which Cottrell-Lomer barriers give way under the shear 
stress due to arrays of dislocations piled up against them. Table 2 
shows this to be the case except at 77°K. (The product y,7'/? shows 
the same trend, indicating a linear relation between og and y,.) Stroh’s 
(1956) more detailed analysis of the strength of Cottrell_Lomer barriers 
shows however that this simple criterion cannot apply at low temper- 
atures. The approximate constancy of o,7'/? would however also be 
expected if transition from stage IT to stage III were due to the circum- 
vention of Cottrell-Lomer barriers by dislocations, facilitated by cross- 
slip (Seeger et al. 1957); the numerical value of the product does not, 
at present, enable one to discriminate between the two processes. 


Table 2 
Xp le 
Of! CED ate Xue ANSE se OND (kg/em2)? 
ok (kg/om?)°KY/? | (kg/em?)?°K1/2 | ke/em?® = |———+________ 
Measured | Eqn. (8) 
77 7-0 5:8 3-6 6-7 6-2 
195 9-4 6-7 3-2 4:8 4-9 
293 9-4 6-2 2-8 3-6 3:8 


The observations of Blewitt et al. (1955) that single crystals of copper 
with orientations favouring slip on a conjugate glide system harden more 
rapidly in stage IT than crystals with orientations away from the [100]- 
[111] great circle seems to show that Cottrell-Lomer barriers are being 
formed also in the process of hardening, and not only at the beginning of 
stage II as originally proposed by Friedel. This point of view is also 
expressed by Diehl (1956) and Seeger et al. (1957). Schmid’s Law is 
thus not strictly true in face-centred crystals, and allowance for this 
must be made in deriving work hardening curves of polycrystals from 
single crystal data on the basis of theories in which Schmid’s Law is 
implicit, e.g. Sachs’s (1928) and Kochendorfer’s (1947). 
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Now, the coefticient of shear hardening in stage II, derivable from the 
curves for copper single crystals of various orientations deformed between 
4 and 300°K (Blewitt et al. 1955) lie in the range+ 1-6-2-4 x 10° ke/em?. 
Sachs’s theory requires the coefficient of linear hardening of the poly- 
crystalline copper to be equal to (2-24)? times the shear hardening co- 
efficient of the single crystal, i.e. at most 1-2 x 104 kg/em2, on using the 
highest value obtained from Blewitt’s data, which applies to a crystal 
with an orientation close to [100]. 

Reference to table 2 shows that y, as measured is 2-3 times greater 
than this estimate. Allowing for the possibility that the factor (2-24)? 
may be an underestimate, it still appears as if the complex system of 
stresses acting on the grains of the polycrystal lead to a somewhat greater 
hardening rate in the individual grains than would be found in similar 
grains subject to uniaxial tension, the complexity of the stress systems 
resulting in slip on several ‘ primary’ slip systems. However, in alumi- 
nium and dilute aluminium alloys Jaoul (1957) found a somewhat 
better correlation between the shear hardening rates of single crystals 
and the hardening of corresponding polycrystals in tension, also on the 
basis of the factor (2-24)?. 

The shear strain over which stage II extends in single crystals is also 
considerably greater than the tensile strain up to og in polycrystals ; 
but o, and the stress 27,, where 7, is the shear stress at the end of stage IT 
in single crystals, are approximately equal (I), as is to be expected from 
the identity of the mechanisms leading to the onset of stage III. 


4.3. Parabolic Hardening 


It was pointed out in I that certain features of the distribution of 
dislocation pile-ups in alpha brass (Jacquet 1954) and the similarity 
between alpha-brass and copper suggested to the present authors that 
hardening in copper in stage III might be of the interaction type which 
leads to a parabolic work hardening law (Mott 1952). The stress-strain 
curves of copper polycrystals in stage III were in fact shown to be sections 
of parabolae, the latter having their origins close to the origin of the 
~ coordinate system, as the parabolae fitted to the curves in the figure. 
Also, an attempt to apply some relations derivable from Mott’s theory 
lead to the conclusion that hardening early in stage III was probably 
due to interaction between pile-ups which had originated from the small 
rectangular slip zones characteristic of hardening in stage II at stresses 
above o,, when Cottrell-Lomer barriers ceased to be effective in confining 
slip to the interior of such zones. The additional experimental results 
obtained in the present work permit a somewhat more detailed examina- 
tion of this hypothesis. In deriving the parabolic hardening law we 
shall adopt the method used by Mott (1952), with a different interpretation 


+ This was erroneously stated to be 1-6-2-4x 10* kg/em? in I. 
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however of the origin of the pile-wps responsible for hardening. We shall 
make the following assumptions : 


(a) At the end of stage II the structure of any grain is as given by 
Friedel’s (1955) model, i.e. a random distribution of rectangular slip 
zones of linear dimensions 2Ly=10, locked on all sides by Cottrell- 
Lomer barriers; on account of their size interaction between them is 
small (Friedel 1955). 

(b) As the stress og is reached the non-uniformity of the internal 
stress will result in the breakdown or circumvention of Cottrell-Lomer 
barriers containing the dislocations of some of these slip zones. 


(c) The pile-up in any given zone which has become active in this 
manner will expand to give a ‘ large’ slip zone of linear dimensions 2L 
and slip distance L. 

(d) The stress relaxation in such a pile-up, resulting from the expansion, 
will facilitate renewed locking of the expanded slip zone by Cottrell— 
Lomer barriers. 

(e) The ‘ large’ slip zones cannot be regarded as dislocation dipoles 
(Friedel 1955), and their long-range stress fields, which oppose the applied 
stress, necessitate an increase in stress above o, if further ‘ large ’ slip 
zones are to be formed. 


(f) The major source of hardening is the increase in the density of 
‘large’ slip zones with increasing stress, at least in the early stages of 
stage III here under consideration ; hardening due either to any increase 
in the slip distance LZ, or in the thickness of these L-zones due to band- 
formation—if such effects occur—is negligible. 


(7) The number of dislocations, n, in the L-zones is constant and equal 
to that in the L,,-zones, such that (Mott 1952) 


N=T(og—o,)Ly/Gb, 7) ae? hs 7a | hace tem eel 


where we have written jo, and 3a; in place of the shear stress tg, and the 
initial yield stress in shear 7,of a single crystal respectively. Gis the shear 
modulus, 6 the Burger’s vector, and the slip distance in an Ly-zone has 
been taken to be equal to Ly. 

We shall also make the simplifying assumption that the internal stress 
pattern in any grain above og is the same as would have been obtained 
if interaction hardening only had taken place right from the beginning, 
i.e. starting with o=0. This assumption amounts to replacing the whole 
hardening curve (e.g. the lowest curve in the figure) by a parabola passing 
through the origin. No physical significance is, of course, attached to 
this model at stresses below og. It also implies that as a small stress 
is applied to the hypothetical model in its undeformed state a 
number of fully formed L-zones will appear instantly, with an average 
spacing between adjacent bands given by x. As the yield stress of the 
undeformed material is assumed to be negligible (the curve passing 
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through the origin) the internal stress o; due to a random distribution 
of L-zones of this type, given by (Mott 1952) 

Garber (4) 
must be exceeded before further hardening can take place. Thus the 
internal stress is equal also to the external applied stress o which must 


be exceeded before additional L-zones are formed. 
The tensile strain corresponding to o is given by 


CU Naeem sd gt yt tO) 


where, as in I, we take the tensile and shear strains to be equal. From 
eqns. (4) and (5) we have 


Gray eC, ae een ay ee Me) OT, fy, 4x0) 
where 
Xp=nGb/n7L, Pee ahi eye eyehet 8. IPT) 
or, on using eqn. (3): 
i= Gi LILe Nope) eA 48) 


The values of y,, evaluated from eqn. (8) with L/Ly =2 and o,=60 kg/cm? 
(Blewitt et al. 1955) are given in table 2 together with the measured 
values. The agreement is satisfactory both from the point of view of 
the order of magnitude, as well as in relation to the temperature-depen- 
dence. (The linear relation between y, and ov, already referred to above, 
is also apparent from eqn. (8).) The coefficient of shear hardening of 
single crystals in stage III should be equal to about (2.24)-% times x,. 
This estimate is in fact of the right order (1). 


4.4, Cottrell-Lomer Barriers 
By using the measured values of og and y, respectively it is possible 
to deduce from eqn. (3) and eqn. (7) two independent values of the 
number, 7, of dislocations which will lead to the break-down or other 
mode of failure of Cottrell-Lomer barriers. Both estimates, together 
with their mean, are given in table 3. 


Table 3 
Gp n 
ite R min max R/ Rin R/Rmax 
Eqn. (3) | Eqn. (7) | Mean 
Wi 88 91 90 0-070 | 0-041 0:20 1-7 0-35 
195 LO Us, 73 0-049 | 0-031 0-15 1-6 0:33 
292 62 61 62 0-035 | 0-027 0-13 1-3 0-27 


Now, Stroh (1956) has plotted the ratio nog/2@ as a function of the 
temperature for the two possible modes of breakdown of a Cottrell-Lomer 
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barrier in copper. The maximum value of this ‘ breakdown ’ ratio 
Roaxs given by Stroh, relates to breakdown by ‘ dissociation =) the 
width of the core of a dislocation is taken to be 0-2b. The minimum 
R, yin relates to breakdown by ‘ recombination ’, the assumed core-width 
being 0:3b. The values of the ratio, denoted by Rin table 3, obtained with 
n and a, as given in table 3 and table 2 respectively, lie between the 
limits of Stroh’s estimates (table 3). The present results do not therefore 
rule out the possibility that a ‘ breakdown’ mechanism of the type 
visualized by Stroh may be operative at the end of stage IT as well as 
in stage IIT. 
However, Seeger et al. (1957) have recently suggested that as 


(a) the activation energy for cross-slip is less than that for ‘ breakdown ” 
and as 


(b) the breakdown mechanism would be expected to result in the 
appearance of long slip lines on the surface of crystals entering stage III, 
a phenomenon not observed by them in copper, dislocations piled up 
against Cottrell-Lomer barriers would surmount the barriers by cross-slip 
rather than break through them. These arguments do not appear to 
us to be fully convincing for the following reasons : 

(a) The relation between the magnitudes of two activation energies 
is not a decisive criterion for discriminating between two processes in 
the absence of a detailed model of the mechanisms of both ; such a model 
is not yet available for the mechanism involving cross-slip. 


(6) Dislocations which have broken through a Cottrell-Lomer barrier 
could either again be stopped by a new barrier, as suggested in this 
paper in dealing with hardening in stage III, or they could be deflected 
to other slip-planes by cross-slip. Cross-slip, occurring subsequent to 
‘break-through ’ could thus account for the fragmentation of slip-bands 
observed for example by Seeger et al. (1957) in copper, and by the present 
authors (Feltham and Meakin 1957 b) in single crystals of lead, while 
the break-through process itself could provide a rational explanation 
of some characteristics of slip in metals, e.g. the occurrence of slip-ledges 
several thousand A high in metals with low stacking-fault energies 
such as alpha-brass (Wilsdorf and Towrie 1956), or the sudden appearance 
of fully-formed slip steps about 2004 high early in the process of band- 
formation in aluminium (Kuhlmann 1950). 

It may well be that within a certain domain of the variables, i.e. 


temperatures, stress and strain rate, one of the processes may be operative, 
outside that domain the other. 
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ABSTRACT 


The specific heats of the metals titanium, zirconium, and hafnium have 
been measured between 1-2 and 20°x. Values of y, the electronic specific 
heat coefficient, have been determined for these metals as 8-5, 7:25 and 
6-3 x 10-4 cal/mole deg, while the corresponding values of the Debye @ in 
the continuum region are 430°, 310° and 261°. The variation of 6 with 
temperature has been investigated, and it is shown that there is a simi- 
larity in this variation for all three elements when a suitably reduced 
scale of temperature is used. 


§ 1. INTRODUCTION 


THE availability of high purity titanium, zirconium, and hafnium in recent 
years has made possible the investigation of the thermal properties of these 
metals. As the information on these elements is not complete below 20°K, 
it was thought of interest to extend the specific heat measurements from 
liquid hydrogen to liquid helium temperatures with the aim in particular 
of determining the electronic specific heat of hafnium which has not 
previously been reported. 

The elements titanium, zirconium, and hafnium in Group IVA of the 
periodic table all crystallize in the hexagonal close-packed structure with 
almost equal axial ratios (c/a=1-59). For close-packed spheres c/a would 
be 1-633. All three elements also have body-centred cubic modifications ~ 
which are stable at high temperatures. Some similarity in their thermal 
properties may therefore be expected. 


§ 2. PREVIOUS MEASUREMENTS 


The specific heats of titanium and zirconium have been measured by 
Johnston and Kothen (1953) and Johnston and Skinner (1951) respectively 
between 15° and 300°x. Other measurements down to liquid nitrogen 
temperatures are those of Kelley (1944) on titanium, and Todd (1950) on 
zirconium. Johnston’s work has been used here for comparison as he 
had the purer sample available, and the measurements extend over a 
wider temperature range. For hafnium only the results of Simon and 
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Cristescu (1934) are available, but it is now thought that their sample 
was largely hafnium oxide. 

At liquid helium temperatures there are only the measurements of 
Estermann e¢ al. (1952) on the specific heat of titanium and zirconium. 
Recently some measurements have appeared on titanium between 4° and 
15°K (Aven e¢ al. 1956). 


§ 3. EXPERIMENTAL DETAILS 


The samples of the elements were provided by A.E.R.E. Harwell and 
were all made by the van Arkel process in which the crystals are formed 
by condensation from the iodide. It is possible that the zirconium 
contained 1°%, hafnium, but the hafnium was said to be zirconium free. 
Other impurities were stated to be less than 0-1°%. The samples consisted 
of about one mole of material for each element. 

The calorimeter used in these experiments has been described in 
detail by Smith and Wolcott (1956). Since the technique of measure- 
ments is also described by these authors it will not be further elaborated. 


§ 4. RESULTS 


In the helium region the specific heat can be described by the formula 
Ch=yTl +4647%/6 cal/mole deg where the linear term in 7’ is the electronic 
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Plot of C,/T against 7? for titanium, zirconium, and hafnium. The ordinate 
is plotted in units of 10-4 cal/mole deg?. 


specific heat and the cubic term the lattice contribution. Above 4° the 
Debye characteristic temperature, #, decreases with rising temperature for 
all three elements. The following values of y and @ were appropriate for 
these substances below 4°x : titanium, y=8-5 x 10-4, 0=430° ; zirconium, 
y=7-25X 10-4, 6=310°; hafnium, y=—6-3 x 10-4 cal/mole deg?, 0= 261". 
Graphs of C,/T vs 7? from which these values are derived appear in 
fig. 1, while the experimental points are given in the table. Since the 
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lattice contribution is less than 10°% of the total specific heat of titanium 
and zirconium in the liquid helium region, @ for these elements can only 
be specified to +3°% below 4°K. 


§ 5. DIscussION 


5.1. Lattice Specific Heats 


In order to determine the true lattice contribution to the specific heat 
and hence the corresponding Debye 8, it is necessary first of all to correct 
the measured values of C,, to C, and then to allow for the electronic term, 
y7. As an indication of the relative magnitude of these corrections, 
neglecting the electronic term for zirconium at 150°K or 7/@ about 3/5 
would introduce an error of 10% to the Debye §, while ignoring the 
C,-C, correction would add another 4%. The correction for the 
electronic term is easily made using the value of y determined in the helium 
region where C=C, and the lattice contribution is proportional to 7°. 
The C,-C, correction can be determined from the approximate relation 
C,-C,=AC,?T where A is determined by comparison at room temper- 
ature with the exact thermodynamic expression C,-C,=«?V7'/8, where 
a is the bulk thermal expansion coefficient, 8 the compressibility, and V 
the atomic volume. For zirconium the value of A was determined to be 
8-9 10-§ mole/cal using the room temperature values of the thermal 
expansion and compressibility (van Arkel 1939, Reynolds 1952). Since 
expansion data are not available for titanium the empirical formula of 
Nernst and Lindemann (1911) was used in which A=0-0214/7',, where 
T, is the melting temperature, giving A=1-07X10-> mole/cal. The 
Nernst-Lindmann formula also gives the value of A=8-9x10-* for 
zirconium, so that the value obtained for titanium in this way is expected 
to be reasonable. Below 20°K, C.-C, is negligible, and so no correction 
has been made for hafnium in this temperature range. 

The 6-T' relations for these three elements are shown in figs. 2 to 4. 

The results are in general accord with the present theory (e.g. Blackman 
1955) from which a temperature variation of the Debye 6 may be expected 
at temperatures above about 6/50 where the non-continuum portion of 
the vibration spectrum becomes reasonably excited. For titanium 6 does 
not change very much with temperature up to 13°K ; this is in agreement 
with the results of Aven et al. who find a value of 6=421 between 4° and 
13° while the present measurements indicate an average 6 value of 425° 
in this temperature range. The present results connect quite well with 
those of Johnston, as figs. 2 and 3 indicate. Aven commented that the 
# values calculated from Johnston’s data seemed much too low ; however, 
he apparently neglected to allow for the electronic specific heat when calcu- 
lating the Debye 6 from Johnston’s measured values of Cy. When correction 
is made for the electronic term the higher values of 0 shown in fig. 2 are 
obtained. On the other hand the @’s in the 2°-4° region given by Ester- 
mann ¢ét al. for titanium (280°) and zirconium (265°) appear to be much 
too low compared with the present value of 430° (Ti) and 310° (Zr). 
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However, these authors also obtain @ values which appear low for copper 
(Corak et al. 1955) and magnesium (Smith 1955), so that the discrepancy 
systematic one. 
ee pe in the 0-7’ curve of these elements is further brought out 
by fig. 5 in which the reduced values 4/ 6,, and T'/6,, have been plotted. 
Here 6., is a characteristic value of 6 at high temperatures. For titanium 
6. was taken as 360°, for zirconium 260°, and for hafnium 200°. (The 
results of Friedberg (private communication) have been used in deter- 
mining 0,, and the shape of the reduced @ curve at high temperatures for 
hafnium.) Exact coincidence in the curves is not to be expected however, 
for as Hill and Smith (1953) have pointed out, this would imply that the 
ratios of the interatomic force constants were identical as well as the 


axial ratios. 
Fig. 5 
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Reduced @ plotted as a function of reduced temperature for Ti, Zr, and Hf. 


The reduced 6 curves may be compared with those for the hexagonal 
metals of Groups II and IIB of the periodic table (Hill and Smith 1953, 
Smith 1955, Smith and Wolcott 1956). When this is done it appears that 
there is a similarity among the Group II elements beryllium and mag- 
nesium, and the Group IVA metals, all of which have axial ratios slightly 
less than the ideal \/8/3=1-633. The temperature variation of @ for 
zinc and cadmium (c/a=1-9) is much more pronounced. One might be 
tempted to conclude from this that for the hexagonal metals large axial 
ratios are associated with large temperature variations of the Debye 0. 
However, rhenium with an axial ratio of 1-6 has a Debye @ which changes 
over 60% below 20° (Wolcott 1955), so it appears that other variables 
as well as the c/a ratio must be significant in determining the @ behaviour 
of the various hexagonal metals. 

It is interesting to note that the high temperature 6’s closely follow the 
Lindemann (1910) formula which requires that @ be proportional to 
(7',/M)'?2V—"3, where 7’, is the melting temperature, M the atomic mass, 
and V the atomic volume. For the elements titanium, zirconium, and 
hafnium this quantity has the proportion 1 ; 0-70 ; 0-54, while the actual 
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high temperature 6’s are in the proportion 1; 0-72: 0:55. This corre- 
spondence might be expected in view of the similar physical properties 
of these elements. 

As no complete theoretical calculations of the vibration spectra have 
been carried out for the hexagonal close-packed lattice, no detailed 
comparison of the @—-T curves with theory can be made. However, Barron 
and Domb (1955) have made some calculations using a central force model 
and conclude that the Debye @ at 0°K shifts about 1°4 when the lattice 
structure changes from hexagonal to cubic close-packed. In view of this 
agreement at low temperatures there may be some justification in compar- 
ing qualitatively the experimental 6-7’ curves for the h.c.p. metals with 
the theoretical calculations for the cubic close-packed structure. When 
this is done, it is seen that the 20°% drop in 6 with rising temperature 
is indeed similar to that found by Leighton (1948) for the f.c.c. lattice. 
It should be mentioned though that there are f.c.c. metals (thorium) as 
well as h.c.p. metals (rhenium) which show much larger @ variations than 
those indicated by Leighton. When the single crystal elastic constants 
for these metals become available a more detailed check of the adequacy 
of the Leighton model can be made. 

A @ may also be calculated from the isotropic Debye model using the 
room temperature longitudinal and transverse wave velocities of Reynolds 
(1952). When this is done values of 6=380° for titanium and 6= 264° for 
zirconium are obtained. The agreement with the high temperature 
@’s (which also occurs in the case of copper) must, of course, be considered 
fortuitous, as the continuum theory is not applicable in this temperature 
region. 


5.2. Electronic Specific Heats 


The electronic specific heats may be related to the density of electronic 
states at the Fermi surface, g(é), since quite generally y=2/377kq(é). 
The present values of y give the density of states to be 0-755 levels/atm ev 
for titanium, 0-644 levels/atm ev for zirconium, and 0-56 levels/atm ev 
for hafnium. These results are about two and one half times the free 
electron values assuming four free electrons per atom, and are therefore 
qualitatively in agreement with Mott’s (1935) theory of the transition 
metals. However, as no theoretical calculation of the density of states 
has been made for any of these elements, no real comparison with theory 
is possible. From the theoretical density of state calculations that 
have been made for other transition metals, the value of y for the elements 
in Group IVA would be expected to be larger than the free electron value, 
but smaller than that for the later transition elements in Groups V and 
VIL, which is in fact what is observed. In common with other transition 
metal series, it appears that the heavier elements have the lower density 
of states; in fact, for Group IVA the density of states is very nearly 
proportional to Z~1\/4, where Z is the atomic number. 

The actual values of the electronic coefficient y obtained here for 
titanium and zirconium are a few per cent higher than those reported by 
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Estermann et al. and Aven et al., but the discrepancies are probably 
within the combined experimental error. 

As titanium and zirconium are superconductors, estimates of y 
may be obtained from their magnetic threshold field curves. For 
titanium the measurements give y=35 x 10-4 (Hein et al. 1957), and for 
zirconium y=4 x 10-4 cal/mole deg? (Smith and Daunt 1952). These 
values are derived from theoretical formulae which treat the critical 
magnetic field curve as representing a phase boundary between the normal 
and superconducting states of the metal. This picture is idealized in 
as much as it admits no intermediate state. In the case of hard super- 
conductors such as titanium and zirconium which have broad transitions 
and therefore large intermediate states, the application of thermodynamics 
to the critical field curves is questionable, and hence the disagreement 
with calorimetric values of y is not surprising. 


§ 6. CONCLUSION 


The low temperature heat capacities of the Group [VA metals have now 
been fairly thoroughly investigated. The electronic specific heats and 
Debye 6’s are in agreement with the general theoretical picture, whilst 
there is a similarity in the thermal properties of these elements which 
is not unexpected in view of their similar physical and chemical properties. 
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ABSTRACT 


The angular correlation of the gamma cascade in 348 has been measured 
to be 1-++ (0-332 + 0-032)P, (cos @)—(0-017-+-0-011)P,,(cos @) and this is inter- 
preted in terms of a spin-sequence 2+, 2+,0 with a mixture of 1-4°% 
of 2 radiation in the first transition. A new angular correlation appara- 
tus and fast coincidence system is briefly described. 


§ 1. INTRODUCTION 


Tue *4Cl nucleus is the lightest nucleus with 4 =4n-+2 to have the level» 
with total isotopic spin zero (7’=0) higher in energy than the T=1 
level, which forms the ground state. Consequently, the beta decay of 
both the isomeric level (¢;,.=33 min) and the ground state (t1;2.=1-5 sec) 
are of particular interest, and the mode of decay is well established 
(see fig. 1). 

Angular correlation measurements for the gamma ray cascade in 48 
can be used to assign spins to the energy levels. Measurements were 
published for this correlation by Handler and Richardson in 1956 whilst 
the present work was in progress, but it was thought that the information 
was of sufficient interest, and the measurements sufficiently complicated 
by the short half life and by the presence of the abundant annihilation 
radiation, to make an independent measurement and analysis of the 
results desirable. 

This paper contains some details of a new angular correlation apparatus 
built to make measurements on isotopes with half lives as short as 10 
minutes. It includes a new and simple fast coincidence arrangement. 
Measurements of the angular correlation for the gamma rays from *4S are 
considered in association with the other data on the decay scheme to 
assign spins to the *48S energy levels. 


§ 2. APPARATUS 


Two scintillation counters were mounted on the movable arms of a 
table so that they could be set at any angle from 45° (when they were 
in contact) up to 180°. Each counter consisted of a 1} in. x 2 in. sodium 
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iodide (Tl) crystal mounted on an E.M.I. Type 6262, 14-stage photo- 
multiplier, with a lead shield 2 in. long defining the half angle of acceptance 
to 10°. The mounts which held the counters to the table were provided 
with fine adjustments so that each counter and shield could be directed 
accurately at the table centre, where the source was situated. 


Fig. 1 
3+ 32.5min T=O 


ies Oe: ERE ial 
1.46 sec. 
(2+) Jo 
(24) J 
(0+) Jy 


The decay scheme of *4Cl to *48. This decay scheme summarizes the results 
of the following authors: Ruby and Richardson (1951), Ticho (1951), 
Aber and Stahelin (1953), Stahelin and Preiswerk (1953), Peaslee (1953) 
Green and Richardson (1954), Handler and Richardson (1956), Bleuler 
and Morinaga (1956). The spin assignments resulting from the present 
measurements are shown on the left in parentheses. 


To form a fast-slow coincidence system (Bell and Petch 1949), two 
output signals were taken from the 14-stage photomultipliers. The 
first, which came from the 9th dynode, was given conventional ampli- 
fication and was used for pulse height selection and spectroscopic purposes. 
The second signal came from the final electrode (collector) and was used 
for the fast coincidence work. With 1-5 kv applied to the photomulti- 
plier the final stages tended to saturate, particularly for large pulses 
resulting in preferential amplification of smaller signals, and reduction 
of all pulses to approximately uniform amplitude. The negative output 
pulse rose to 0-3 v in 10-8 sec when fed directly into a load of 100 ohms 
and a 110 0hm cable. At the far end the cable was unterminated and 
fed pulses into an E.F.P. 60 secondary emission pentode with an inductive 
load in the anode, which in turn drove a low level trigger circuit using 
another E.F.P.60 pentode. This trigger circuit has been described 
by Moody and McLusky (1952). It has accurate and low level triggering 
because the valve is passing appreciable cathode current in its quiescent 
state and consequently works under approximately optimum conditions 
The feedback from secondary cathode to grid is a direct capacitative 
coupling involving only small strays, and the whole system is kept 
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stable by interrupting the feedback path with a germanium diode just 
biased-off. The circuit was originally designed as a discriminator, for 
which the very low triggering voltage of 0-05 volts gives accurate selection 
of small pulses (0 10v). The same factors which make it suitable 
as a low level discriminator (i.e. sensitive triggering and production of a 
standard output pulse) also make this circuit suitable for use in a coin- 
cidence arrangement. The circuit is triggered by very small pulses 


SINGLE CHANNEL 
KICKSORTER 


COUNTER I 


SCALER. 


COUNTER IT 


to slow amp « 
+ 500v discriminator + 500Vv 


to bias potential 
of O.lv 
Signal from 
fast amp II 


to bias potential 
iy a es 
Signal from Se 


fast amp. | 


TRIGGER CIRCUIT L COINCIDENCE DIODE TRIGGER CIRCUIT DT 


The fast-slow coincidence arrangement showing a detail of the fast coincidence 
system. 


(about 0-05 v) and if the control grid is driven through even a few tenths 
of a volt in 10~® sec the positive feedback, which has a very small lag, 
rapidly drives the valve into saturation current. The triggering point 
of the circuit depends upon the bias across the diode in the feedback 
path and when this bias is reduced to the minimum consistent with 


4Q2 
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long term stability of the circuit, triggering by the signal from the photo- 
multiplier (after amplification in the secondary emission pentode) is 
estimated to be accurate to within 2x 10-®sec. When triggered the 
circuit gives a 5 volt pulse across a 110 ohm line 5 mysec long, shorted 
at the far end to clip the signals to 10~* sec. 

Clearly, two such signals mixed in a coincidence circuit will give a 
resolving time of about 10-8 sec. In the present apparatus, pulses 
from the trigger circuits are fed to a biased-off germanium diode followed 
by slow amplification. The schematic arrangement is shown in fig. 2; 
apart from the details of the fast coincidence circuit, it is a standard 
fast-slow coincidence assembly. Pulses selected for pulse height are 
put into slow coincidence (27=10~* sec); the resulting coincidences 
are formed only between pulses of the required amplitude (and hence 
correspond to the energy range of interest) but a large proportion of 
random coincidences may be included at fast counting rates (10° per sec). 
The fast coincidence system gives coincidences regardless of pulse height, 
but includes a very low proportion of randoms. Coincidences common 
to both systems are selected in a further coincidence circuit, and these 
so called fast-slow coincidences occur only between pulses of the desired 
amplitude and contain the low proportion of random characteristic of 
the fast coincidence circuit. 

The apparatus was checked as a whole using °°Co for which the angular 
correlation is well established and the results of this measurement are 
given in fig. 3 together with the theoretical values; the agreement is 
satisfactory. 


§ 3. OBSERVATIONS 
34C] may conveniently be formed by the reaction 
328 +4He > 4™Cl+ pt+n 
and this reaction was induced by bombarding a natural sulphur target 
with the external beam of 40 Mev He?* ions from the 60in. Nuffield 
cyclotron of this laboratory. The external beam is free from contamina- 
tion by the low energy helium ions and deuterium ions which might be 
present in circulating beams within the cyclotron tank. The sulphur 
target, which was about 2mm thick, adhered to a water cooled copper 
backing and was covered during bombardment by a 0-0003 in. tantalum 
foil. Bombardments of 0-5 ja for 10 minutes provided sufficient activity 
to make up sources of convenient strength. Preliminary experiments 
confirmed the prediction that under these conditions there would be no 
competing nuclear reactions induced which would produce activities 
comparable in intensity with that of the ##"Cl after the lapse of a few 
minutes. Consequently it was possible to use the bombarded sulphur 
directly as the source, without chemical extraction of the 24™Cl activity. 
Sulphur from the target was scraped off and sealed into a thin walled 
(0-002 in.) brass cylinder 2mm in diameter and 5mm long, extending 
from a 2 B.A. brass threaded rod at the lower end. This sealed source 
could then be screwed into the centre of the angular correlation table. 
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Direct use of the target sulphur in the source holder had the great advan- 
tage of shortening the time between the bombardement and the measure- 
ments, and also avoided the use of silver chloride from the standard 
chemical extraction of chlorine, with its attendant high density, as the 
source. Initial source strengths used for the angular correlation were 


Fig. 3 


N60 


The angular correlation for Ni. The broken line represents the least squares 
fit. The solid line is the theoretical correlation for a 4(2)2(2) 0 transition. 


of the order, of 0-1 mc and these gave a ratio of real to random coin- 
cidences of about 20 with the resolving time of 10-%sec. The pulse 
height analysers of both scintillation counters were set to receive pulses 
corresponding to an energy dissipation in the crystal in the range 0-7 to 
2-5 Mev. 

The half life of 34"Clis only 33 min, and to make the angular correlation 
measurements it was necessary to replace the source every 80 minutes 
and then to relate the measurements together. It was decided to take 
measurements with each source at 15° intervals for 7 angles for a period 
of 10 minutes at each angle. At each observation a period of 1 minute 
had to be allowed for reading and re-setting the scalers. For a given 
number of observations in a run this reading time was fixed, and to 
concentrate it in early stages of the measurements when the source was 
most active, which would have been the result of taking measurements 
after a constant count had been recorded, meant losing a great deal 
more information than in counting for constant time. The total number 
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of coincidences at angles with cos?@ less than 0-40 was made roughly 
equal to the total number at angles with cos*@ greater than 0-40 by 
taking observations in the following sequence : 
90°, 165°, 105° or 75°, 150°, 120° or 60°, 135° or 45°, 180°. 

With each new source a start was made at the angle next in the sequence 
to the one which started the previous series. 

Frequent checks were made during the course of the experiment on the 
resolving times of the coincidence circuits. 


§ 4, TREATMENT OF RESULTS 


The observed counting rates had to be corrected for the dead times 
of the circuits, which were determined in all cases by the 5 usec dead 
time of the 1009 scaling units, and allowance had to be made for the 
inclusion of random coincidences in the coincidence count. ‘To compen- 
sate for any slight decentring of the source, and small variations in the 
discrimination levels, the combination Y=(C/N,N,.) exp (—At) was 
formed. CO was the coincidence rate after subtraction of the calculated 
random coincidences, N, and N, were the single counting rates in the 
two counters after correction for the scaler dead time, and A was the 
decay constant of Cl. The quantity Y is not affected by the decay 
of the source, although it is strongly dependent upon the source strength 
at time t=0, and consequently relative values of the function at different 
angles are directly proportional to the correlation function W(@) for 
those angles. Observations from different sources were related together 
by fitting the observations to a two parameter curve Y=a)+a, cos? @ 
and then using the a) as a normalizing factor. The results of the observa- 
tions could then be collected together, averaged and fitted to a three 
parameter curve. 

ay+a, cos? 0+-a, cos* 6 


by the method of least squares. The observations at 180° were not 
used in the evaluation of the normalizing factor, or when fitting the 
observations to the curve, since the presence of annihilation radiation 
was likely to induce at this angle a rise in the coincidence rate which 
could not be readily estimated. 

The (x?) significance test (Kendall 1945) was applied to the deviation of 
the observations from the least squares fit taking care to make allowance - 
for the various constraints imposed during the calculations. The test 
gave a value for the chance deviation from the calculated curve as 80%, 
indicating that all the parameters ao, a, and a, were significant. Having 
thus established the validity of the parameters, the errors in the co- 
efficients a), 42 and a, were computed from the mean square deviation 
of the observations from the fitted curve. 

Kach counter had a small but significant efficiency for the exceedingly 
abundant radiation resulting from the annihilation of the positrons. 
Much of this radiation was coincident with the various gamma rays of 
the decay, and it was necessary to correct the observed distribution for the 
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inclusion of the isotropic contribution from gamma-annihilation coin- 
cidences. This contribution was estimated by carrying out a correlation 
experiment at 5° intervals from 150° through 180° to 150°, with each 
spectrometer in turn set upon the annihilation radiation photopeak. 
From the rise in the coincidence rate at 180° it was possible to estimate 
the relative efficiency of the unaltered counter for annihilation radiation. 
In this way the annihilation gamma contribution to the isotropic coin- 
cidence rate was estimated as 8-1-.0-5%. 


Fig. 4 


60 90 120 . 150 180 
; e— 
The Angular correlation for *48. The curve drawn through the observations 
is the least squares fit. 


Finally, the result was converted to the standard form : 
W(0)=A,+A.P,(cos 6)+A ,P,(cos 6) eee) 
and the coefficients A, and A, corrected for the finite solid angle sub- 
tended by the detectors. This correction has been treated by Rose (1953), 
and his correction factors Q for the geometry used are Q=1, Q.=0-97, 
and Q,=0-88. The experimental results after correction by all these 
factors then appear in the form (1) with 
A,=1, A,=0-3320-032, A,=—0-017+40-011. 
These results may be compared with those of Handler and Richardson 
(1956) if the latter are expressed in the same form (1); they give 
A,=1, A,=0-330+40-055, A,=—0-018+0-021. 


§ 5. Discussion 


The presence of a strong crossover transition in the decay scheme 
of the even-even nucleus *48, together with the knowledge that the 
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3-32 mev: level does not have a long lifetime (<10~7 sec), limits the 
spins which need be considered for the two excited states to 1, 2 or 3. 

The experimental coefficients have been compared with the predictions 
from all permutations of these spins using the tables of Biedenharn and 
Rose (1953). Since a mixture of multipoles is possible in the first transi- 
tion the mixing parameter § has been introduced in accordance with the 
definition of Biedenharn and Rose. For a mixed transition the correla- 
tion function W(#) may be written : 


W(0)= W (0) +0? W ,,(8)+ 26 W 4 11(8) 


where W,(@) and W,,(9) are correlation functions for the pure cascade, 
W41,(0) takes into account interference between W,(6) and W,,(@), and 
82 is the mixing ratio. Values of 5 have been found which bring the 
theoretical coefficient A, of eqn. (1) to the experimental value (A ,=—0-332). 
In general, two solutions for 5 are obtained for each spin sequence, and 
to each 8 there is a corresponding value for the coefficient A,. In the 
table the spin sequences considered are listed with the corresponding 
values for 6 and A,. It is interesting to note that even with mixing 
the spin sequence 3, 2, 0 cannot be reconciled with the observed value 
for A,. If the theoretical coefficient A, is required to le within two 
standard deviations of the experimental values, only those entries in 
italic type in the table are acceptable. 


The Values Tabulated for the Mixing Parameter 6 give A,—0-332 for 
each Spin Sequence. The Predicted Value of A, is shown beside 
the Corresponding Value for 6. 


1; J, ohn 5) A, ro Ars 

] 1 0 —0-477 0) —2-10 0 

2 1 0 0-445 ®) 7:67 0 

3 i 0 0-72 0 op 0) 

] 2 0 0-564. —(0-184 6-73 —0-745 

2 » 0 0-12 0-0046 1-67 _ 0-240 

3 2 0 no real solutions 

1 0 0-122 —0-0407 3:37 0-0764. 

: 5 0 0-411 0-0207 6-90 0-140 
@) 


—0-055 —0-00029 1-48 —0-0654 


Experimental values A)>=1, 43;=0-332+0-032, 4,=—0-017+0-011. 


The assignments with spin 1 for the first excited state all require 
the postulate of strong mixing and it is hard to see how assignments 
1,1,0; 2,1, 0 or 3, 1,0 could be reconciled with the observed mode of 
decay from the isomeric level in *4Cl for which the 3+ assignment is well 
substantiated (see fig. 1). The assignments 3, 3,0 and 2, 2,0 remain. 
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For octopole crossover radiation (y, of fig. 1) to compete with the 
dipole transition y, in the spin sequence 3, 3, 0, it must be E3 competing 
with M1 in a sequence 3—(—=odd parity), 3—,0-++. On this assign- 
ment beta decay from the 3+ isomeric state of °4Cl to the first two excited 
states of *4S would be first forbidden, whereas the measured log ft values. 
indicate that they are allowed. This spin sequence may therefore be 
excluded. 

The spin sequence 2, 2, 0 requires the assumption of only slight mixing 
(0? =1-4%,) to bring the experimental and theoretical angular correlations 
into agreement. With a sequence 2+,2+,0-+ the mixture would be 
M1 with E2 radiation in the 1-2 Mev transition between the levels of 
spin 2, followed by an E2 transition to the ground state. The crossover 
radiation from the 3-2 Mev 2+ level would also be E2 and this would be 
expected to compete favourably with the M1+E2 mixture. Beta 
decay to the 2+ levels from the 3+ isomeric state of #4Cl would be 
allowed, and this agrees with the observed log ft values. The assign- 
ment 2-+ for the first excited state is expected from comparison with 
similar nuclei and although mixing has to be introduced in order to fit 
the observed angular correlation, there appears to be no ambiguity 
in interpreting the results. The experimental and theoretical values for 
the coefficient A, can only just be reconciled however, and it is hoped 
soon to confirm this assignment by making polarization correlation 
measurements with the polarimeter now in use in this laboratory. 
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SUMMARY 


The etch-patterns due to helical dislocations are discussed. It is 
suggested that. under certain conditions the etching of a helical disloca- 
tion gives rise to a spiral etch-pit of the kind observed on germanium 
and silicon crystals. Etch-patterns on (111) faces of fluorite crystals 
containing helical dislocations can readily be interpreted on the same 
basis. 


§ 1. INTRODUCTION 


Ir has been shown previously (Amelinckx, Bontinck, Dekeyser and 
Seitz 1957) that some dislocations can climb into a helical shape. Such 
dislocations have been made accessible to direct observation by various 
decoration processes in crystals of CaF, (Bontinck and Amelinckx 1957, 
Bontinck 1957), NaCl (Amelinckx, Bontinck and Maenhout-Van der 
Vorst 1957, Barber et al. 1957) and TiO, (Van Keymeulen and Dekeyser 
1957). 

Double rows of etch-pits which, on continued etching approach one 
another and finally coalesce, have been described by Bontinck (1957) 
on CaF,. They have been interpreted as etch-patterns of helical disloca- 
tions having their axis parallel to the plane of observation. Other 
patterns which we interpret as being due to the etching of helical dis- 
locations having their axis oblique with respect to the surface, have 
since been obtained. It is the object of this paper to discuss their general 
aspect for any orientation with respect to the observation plane. 


§ 2. THEORETICAL 


We consider first the ideal situation whereby the axis of the helix 
is perpendicular to the surface. During the etching process the pit 
will at every moment be centred on the dislocation line, and penetrate 
in depth at a constant speed v,. Viewed from above, this means that 
the centre of the etch-pit describes a circle during the etching. Sideways 
etching proceeds simultaneously and uniformly at a speed v,, once a 
pit is started. The resulting shape of this combined dissolution can 
roughly be described as a ‘helical cone’. Formulated more precisely, 
one can say that the surface of the pit is generated by a circle, lying 
in a plane parallel to the plane of observation, whose radius decreases 


y+ Communicated by the Authors. 
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linearly with depth and whose centre describes the helix, Figure 1 
represents the projection of such surfaces on a plane parallel to the axis 
of the helix, for three different ratios of v,/v,. 


Fig. 1 
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Projection of a spiral etch-pit on a plane passing through the axis of the helix, 
which is normal to the plane of observation. 


It is possible to formulate quantitatively the characteristics of the 
resulting feature. If the dislocation line is given by 


ef s(t) Mn tae rey Mek +1) 
y=f(t) ee ee aes ie hee FC) 
Cee eee ee ee. (8) 


the equation for the surface of the etch-pit,; assuming the observation 
plane to be (xy), can be obtained by eliminating ¢ between (3) and (4) 


[e—fi(P+[y—folt) P= (iene. eee 4) 


In the specific case of the helical dislocations we can take as a parameter 
equation for the dislocation line 


g==Rcososs y=Rsing; 2=—pd/2r . . - - (5) 
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@ being an azimuthal parameter ; & the radius of the helix and p the 
step. It is easy to see that the projection of the etch-pit surface on the 


plane wz is 
2=R cos (27z/p)(v,/0n)e-- - - « + + (6) 


This is the curve shown in fig. 1. The equation for a section with a 
plane through the axis of the helix is 


a=R cos (272z/p)+[(v,/v,,)?22—R? sin? (27z/p) 7... (7) 
For z> (2) and (3) become equivalent. 
Fig. 2 


Shape of the monolayer step attached to the screw component of the emerging 
poe dislocation. Bunching of steps gives rise to a macroscopic 
spiral. 


The lines seen when observing from above are evidently due to the 
steps marked by shading on fig. 1. We assume that the precise location 
of the step is in the points P;. The spiral as observed (fig. 1) is then 
evidently the geometrical locus of the projection of all such points, for 
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all possible cross sections. A simple, straightforward reasoning, shows 
that its equation in polar coordinates p and ¢ is: 

p=R—ad a=pv,/2r0,,. 
The spiral has Archimedian character, the distance d between successive 
turns being a constant : d=pv,/v,. The step height of the spiral turns 
out to be simply the step of the helix and it has, as a consequence, no 
direct relation with the Burgers vector of the underlying dislocation. 


Fig. 3 


Polygonized spiral etch-pit. 


From fig. 1 it is evident that the step, defining the spiral, becomes less 
and less well defined as v,/v,, increases ; one should observe a decreased 
contrast for such cases. 

For small values of v,/v,, (curve 1) the cross section presents concave 
parts, and the centre of the spiral will not be visible. For large values 
of v,/v, (curve 3) the steps which define the spiral are ‘smeared out ’. 
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The conditions for the occurrence of a vertical segment in the step is 
roughly pv,<27Rv,. The centre of the spiral now becomes visible. 
For intermediate cages one can have rather well defined steps perpendi- 
cular to the plane of observation (curve 2). 

When focusing attention on the unwinding of the monolayer step 
attached to the screw component of the dislocation, one finds that this 
step acquires the shape shown in fig. 2, i.e. a spiral of which the instan- 
taneous centre describes a circle. This causes bunching of monolayers 
in the way shown and the macroscopic spiral now becomes visible. 

If the anisotropy in v, is pronounced one obtains polygonized spiral 
pits of the kind shown in fig. 3, for the case of tetragonal symmetry. 
Polygonization can also occur as a consequence of the polygonization 
of the helix itself (Bontinck 1957, Amelinckx, Bontinck and Maenhout- 
Van der Vorst 1957). 

As in general helical dislocations have their axis parallel to a possible 
Burgers vector it is to be expected that the best developed spirals should 
be observed on planes perpendicular to a possible Burgers vector. 

Etching of helices of which the axis is oblique with respect to the 
surface gives rise to patterns which can be deduced using similar reasoning 
as above. The step height is now the projection of p on the normal 
to the surface. The spiral will be deformed by bunching of the macro- 
scopic steps, in the direction which is the projection of the axis of the 
helix on the observation plane. 


§ 3, EXPERIMENTAL 


Octahedral cleavage faces of CaF’, crystals of which it was established 
previously by means of decoration methods, that they contained helical 
dislocations, were etched with concentrated sulphuric acid. Although 
the circumstances for the observation of spiral etch-pits due to helices 
are not favourable, because of the inclination of the axis of the helix 
(which is parallel to a [110] direction) with respect to the (111) plane, some 
features were observed, which are readily interpreted as being due to 
etching along helices. Figure 5, Pl. 48 shows an example; the spiral 
character is not apparent as a consequence of additional bunching 
of the large steps, in one sector. This is readily explained when con- 
sidering the geometry in detail. The axis of the helix being, e.g. [110], 
the additional bunching, if due to the obliqueness of the axis of the helix, 
has to take place in the direction which is the projection of [110] on 
(111), ie. in the [121] direction. This is precisely what is observed. 

Additional evidence that the described feature (fig. 4) originates from 
helical dislocations was obtained by the study of specimens from crystals 
presenting varying densities of helical dislocations, as deduced by the 
pecoration method. It was found that the number of pits of the shape 
of fig. 5 is proportional to the number of helices detected by the other 
method. It is further evident that the asymmetry of the terraced 
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pits cannot be due to the geometry of the (111) face of fluorite ; this 
would require a threefold axis of symmetry as presented by the majority 
of pits. The sector of bunched steps is moreover found with equal 
probability in any of the three sectors, which is again in agreement with 
the statistical distribution over all [110] directions, for the axis of the 
helices. A direct correlation of a decorated helix and the etch-pit 
centred on it proved to be difficult. 

Attempts to study polished (110) faces were unsuccessful because of 


the difficulty of obtaining good etch-pits on this face. 
[on [07] 
[170] 


Etch-step on a (111) plane of fluorite due to a sequence of prismatic loops. 
The radial discontinuity in the [112] direction results from the particular 
geometry of the ‘twin’ pits. P, and P, are the emergence points 
of the loop. 


Fig. 4 


[112]: projection of [110] on (111) 


Occasionally one observes pits of the shape shown in fig. 6, Pl. 48, 
which present a radial discontinuity. It is thought that these originate 
in sequences of prismatic loops of the kind observed in CaF, (Bontinck 
1957) and in NaCl (Barber et al. 1957) by means of decoration methods. 
It is clear that only sequences having their axis inclined with respect 
to the surface can be etched, and only if the following condition is satisfied 
pl2R<tgx, where « is the angle between the axis of the sequence of 
loops and the normal to the surface. This condition is generally satisfied 
for the CaF, helices. The resulting feature is now a pit consisting of 
closed terraces ; the centre of successive terraces being shifted along a 
[112] direction, bunching of steps will again take place in this direction. 
An additional peculiarity results from the fact that each terrace is 
generated from two centres lying on a line in a [110] direction. As 
etch-pits develop [110] edges this means that the geometry of these 
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‘twin’ pits, is as shown in fig. 4. This explains the occurrence of the 
radial discontinuity in a [112] direction. 

During a recent discussion at Bell Telephone Laboratories W. G. Pfann 
raised the question of the origin of the spirally terraced etch-pits, which 
were first observed in germanium by Ellis (1955), and which have since 
then been observed in germanium and silicon by a number of other 
workers. One of us suggested that they might be due to helical 
dislocations. Although Mr. Pfann’s suggestion has in part stimulated us 
to do the present work, we only have evidence for calcium fluoride. 

The spirals on Ellis’s photographs are very approximately Archimedean ; 
he attributed them to large screw dislocations, the step height of the 
spiral being equal to the Burgers vector of the dislocation. Another 
explanation not based on the presence of dislocations was given by 
Lang (1957). There is however no reason why helical dislocations should 
not be present in germanium and silicon, which are also relatively rapidly 
-cooled melt grown crystals. According to Apker (private communica- 
tion), Dash (1956) using his copper decoration technique, has found some 
indication for helical dislocation in silicon. 
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ABSTRACT 


The tensile yield and fracture strengths of a structural steel have been 
measured before and after irradiation with neutrons. The fracture 
strength increases with irradiation at the same rate as the yield strength 
and the embrittlement cannot be ascribed to a loss of cohesive strength. 
Measurements have been made of the effects of grain size and testing 
temperature on the yield and fracture strengths; applied to Petch’s 
equation they show that the radiation hardening is caused essentially 
by an increase in ‘ frictional ’ resistance in the slip planes, such as would 
be produced by dispersed islands of radiation damage, and that the 
cohesive strength of the lattice and the dislocation locking stress are not 
detectably changed by the irradiation. Notched-impact tests have also 
been made to show the embrittling effect of the irradiation. 


§ 1. IvTRODUCTION 


THE hardening and embrittlement of structural steel by irradiation 
with neutrons (Kunz and Holden 1954, Sutton and Leeser 1954, Wilson 
and Billington 1955) raise several questions. For example, is the em- 
brittlement caused by the hardening, or is there a loss of cohesive strength ? 
To what extent is the hardening caused by the segregation of point 
defects to dislocations, or by the segregation of these defects into clusters 
dispersed through the lattice ? 

We have approached these questions experimentally by using a method 
developed originally by Petch (1953) for analysing the effect of grain size 
on the (tensile) yield and fracture stresses of steel. These stresses, which 
are denoted here by oc, and o, respectively, have been shown by Hall 
(1951) and Petch (1953, 1954) to vary with the grain size / of the metal as 


o,=o,+kl- i, o,=o,+hJ-. Se en ae (1) 


+ Communicated by the Authors. 
t Work done while on attachment at A.E.R.E, 
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Here k,, measures the stress needed to start plastic yielding in the grains 
ahead of the plastic front and varies sensitively with the temperature and 
rate of straining, whereas k, measures the stress needed to crack the 
metal at the tip of a slip band and is almost independent of the tempera- 
ture and rate of straining. The stress o,; measures the ‘ friction ’—due 
for example to precipitates, other dislocations, or the lattice itself—that 
opposes the passage of a free dislocation across a slip plane. The 
theoretical interpretation of these formula is that, when a slip plane is 
made active by an applied stress c, a part of this stress is used to over- 
come friction in the plane, and the rest of it, o—o,, is free to pile 
dislocations against the grain boundary at the end of the plane. The 
stress concentration set up by these dislocations can become large 
enough, when o—<o, is large, to trigger slip planes into activity in the 
next grain, or to open a crack at the end of the slip band. Equations (1) 
give o, and o,, the two values of o at which these effects take place. 

It follows that o,, k, and k,, can be determined from measurements of 
the yield and fracture stresses on specimens of different grain sizes. If 
o, is increased by irradiation we may conclude that the frictional resistance 
in the slip planes is increased by the presence of point defects dispersed, 
as individuals or clusters, throughout the lattice ; and if /, is increased 
then the locking of dislocations is increased by the irradiation. Similarly, 
if k, is decreased then the cohesive strength is decreased by the irradiation. 
In this note we describe an experiment aimed to ascertain these effects. 


§ 2. METHOD 

Specimens were made from bars of En2 steel (0-15% C, 0-6°% Mn, 
0-12° Si) similar to that used by Petch. These were heat-treated to 
give two grain sizes, coarse (500 grains mm*) and fine (2000 grains mm~2), 
and irradiated in aluminium cans filled with helium in the BEpo and 
Windscale reactors to integrated fluxes of 2:5 101" to 5-7 1019 thermal 
neutrons cm? at temperatures of 50 to 100°c. The fast neutron flux 
above 1 Mev was about 75°, of the thermal flux. 

Notched impact specimens were used to show the embrittling effect of 
the irradiation and were tested in a balanced Hounsfield impact machine 
at various temperatures between —60°c and +180°c. The stresses oF 
and o, were measured on standard Hounsfield tensile specimens in a hard- 
beam machine at both room and liquid nitrogen temperatures. 


§ 3. RESULTS 


The figure shows the embrittlement. Small irradiations are able to 
reduce the energy of ductile fracture and increase the temperature of 
brittle fracture appreciably in this steel. 

The table summarizes the tensile measurements. Some specimens, 
particularly the fine grained ones, were ductile at the temperature of 
liquid nitrogen. It is known from the work of McAdam et al. (1945, 
1947), Hollomon and Zener (1944), and Bridgman (1944), that plastic 
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deformation raises the fracture stress. To obtain an undistorted measure 
of the direct effect of irradiation on the fracture stress it is thus necessary 
to separate the indirect effect which the irradiation exerts on the fracture 
stress through the change it produces in the ductility of the metal. We 
have, therefore, used the method of Petch (1953) to reduce the fracture 
stresses given in the table to zero elongation. 


UNIRRADIATED 
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TEMPERATURE °C. 


Effect of irradiation on the impact transition of En2 steel. 
, (Grain size 2000 grains/mm2.) 


§ 4. Discussion AND CONCLUSIONS 


The main feature of the results in the table can be seen by comparing 
the unirradiated and heavily irradiated specimens. Irrespective of grain 
size, temperature, or whether yield or fracture is considered, the strength 
of the metal is increased by about 32 000 lbin~® by the irradiation. 
The only term in eqns. (1) that could give such behaviour is the friction 
stress, o;, and so we conclude that the irradiation has increased this term 
by about 32 000 lb in-?. This strongly suggests that radiation hardening 
in these specimens is caused by dispersed regions of radiation damage, 
such as clusters of point defects, that are sufficiently large to withstand 
thermally-activated movements of dislocations. 

It is clear that the irradiation has not significantly changed the 
cohesive strength of the metal, since the effect of grain size on the fracture 
stress is practically unaltered by the irradiation. The value of k,, about 
2-5 10% c.g.s. units, obtained by fitting the results to eqn. (1) agrees 
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reasonably with the theoretical value 2x 108, deduced by Stroh (1954, 
1955) from the cohesive properties of iron, and also with Petch’s experi- 
mental value. 


Effects of irradiation on the tensile strength of En2 steel. Each value of 
strength is averaged from three or more specimens 


Lower yield stress o, at Lower yield stress o, at 
ena room temperature liq. N temperature 
neutron (Ib in.-2) (Ib in.—2) 
Flux 
os Coarse grain | Fine grain | Coarse grain | Fine grain 
Unirradiated 28 500 36 500 106 000 120 000 
Toe Not tested 123 000 
pn As 
1-5 x 101° 48500 | 54500 143 000 
5-7 x 10 60.500 | 68.000 Fractured ——~ 
Integrated Sea a pay otic °¢ | Reduction in Aréa, per cent 
Neutron pater Ib een ae at liq. N temperature 
Flux Op at a 
=9 
Oey Coarse grain | Fine grain | Coarse grain | Fine grain 
Unirradiated 102 500 124 000 8 42 
Los 115 000 124 000 0-5 36 
Ib 10% 127 000 146 000 0 0-35 
5-7 x 109 133 000 157 000 0 0-20 


It also appears that the irradiation has not significantly changed the 
dislocation locking stress, i.e. the stress to start a slip band, since there is 
practically no change in the effect of grain size on the yield stress. How- 
ever, in the tests at room temperature k, is so small, about 9x 107 ¢.g.s. 
that all but large changes in this term would remain undetected and the 
most we can infer from the present results in that k, remains constant to 
within 15%. We hope to make further measurements at low tempera- 
tures, where 4, is large, to determine whether smaller changes occur. 

It seems reasonable to expect some change in k, in view of the well 
known effect of irradiation in suppressing internal friction due to dis- 
locations in metals such as copper, which has been attributed to the 
pinning of dislocations by point defects (Thompson and Holmes 1956). 
However, the irradiations used in the present work are two to four orders 
of magnitude greater than those at which the pinning effect appears, and 
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under these conditions the change in o, is predominant. It may be 
difficult to detect the effect of pinning on k, in iron, even after light 
irradiations, in view of the strong pinning already produced by carbon 
and nitrogen atoms segregated to dislocations. 
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ABSTRACT 


Crystals of silver halides which have been grown and annealed in air 
have always been used for measurements of photoconductivity. Such 
crystals, which are dilute solutions of silver oxide in the silver halide, 
have a much greater photoconductive response at room temperature than 
crystals of the highest purity. The mechanism responsible for the 
sensitization of the photoconductivity by silver oxide is discussed. The 
relationship of the model which is proposed to that developed by Newman, 
Tyler and Woodbury for the interpretation of sensitized photoconductivity 
in n-type germanium is indicated. 


Tue photoconductive sensitivity of crystals of silver halides is increased 
by annealing them at a temperature near 400°C in air or in oxygen, and 
decreased by annealing them in halogen or in vacuum. The purpose 
of this note is to suggest an interpretation of this behaviour, aspects of 
which have been recognized by all who have studied electronic conductivity 
in silver halide crystals (Toy and Harrison 1930, Hecht 1932, Lehfeldt 
1935, van Heerden 1945, 1950, Haynes 1948, Haynes and Shockley 1948, 
1951, Allemand and Rossel 1954, Brown 1955). 

The experimental basis of the model which will be proposed is provided 
by observations made at Bristol on the displacement in electric fields 
of the photoelectrons which are liberated when light is absorbed in silver 
halide crystals at room temperature. The techniques used were intro- 
duced by Haynes and Shockley (1948, 1951) who mounted air annealed 
silver chloride crystals (Haynes 1948) between two glass plates, the outer 
surfaces of which were coated with platinum to form the two plates of a 
condenser. The front surface of the crystal was illuminated with short 
duration light flashes through gaps in the platinum coating and, simul- 
taneously, voltage pulses appeared across the condenser. The lines of 
force of the resulting electric field were normal to the surfaces of the 
erystal, and the field was in the direction which would cause electrons 
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to drift from the front to the rear surface. After many hours of repetitive 
flashing and pulsing a deposit of silver became visible on the rear surface. 
From this observation, Haynes and Shockley concluded that electrons 
were liberated by the absorption of the light flashes, displaced through 
the crystal by the voltage pulses, and trapped at the rear surface. 

At an early stage in the experimental work on photographic sensitivity 
at Bristol, this method was employed to displace latent images, which are 
formed by much shorter exposures, to the rear surfaces of silver halide 
crystals. The thin single crystal silver halide sheets which were being 
used at that time (Keith and Mitchell 1953 a, b) were probably not very 
different in composition from the air annealed material used by Haynes 
and Shockley (1951), and a high surface density of specks of silver (fog 
specks in photographic terminology) always appeared when they were 
treated with photographic developers without previous exposure to 
light (compare, Haynes 1948, Allemand and Rossel 1954). With such 
sheets, Hedges and Mitchell (1953) were able to produce a developable 
latent image in the areas between the fog specks on the rear surfaces of 
the crystals after about ten light flashes with synchronized applied voltage 
pulses. 

The material used for these experiments was, however, known to be a 
dilute solid solution of silver oxide in the silver halide and the experi- 
mental work was next concentrated on the purification of the silver halide 
(see Clark and Mitchell 1956 for a detailed account). The density of fog 
specks was reduced to a negligible value by methods of purification and 
crystallization designed to eliminate silver and silver oxide from the 
fused silver halide and from the crystals grown from it. The progress 
made can be assessed by a comparison of the surface density of fog specks 
in the photomicrographs reproduced in the papers by Keith and Mitchell 
(1953 a, b) with that in the photomicrographs in the paper by Evans 
et al. (1955). The quality of the large single crystals and thin single 
crystal sheets which were grown from this material, as judged from 
insensitivity to light and the absence of spontaneous reduction on immer- 
sion in active photographic developers, was probably as great as that 
of any previously prepared. 

Evans (1955) then spent many months in attempts to displace the 
latent image at room temperature in crystals grown from this highly 
purified silver halide by the application of electric fields as in the previous 
work of Hedges and Mitchell (1953). He used precisely the same 
apparatus for producing short duration light flashes and synchronized 
voltage pulses but was finally forced to conclude that it was not possible 
to displace photoelectrons at room temperature through well annealed 
crystals of pure silver bromide, free from silver and silver oxide, to 
produce a developable latent image on the rear surface. A large number 
of these experiments were carried out with sheet crystals, 150 in thick- 
ness which had been annealed, first in vacuum to break down any residual 
silver oxide, then in bromine to convert any silver present to silver 
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bromide, and finally in vacuum again to remove excess bromine. If 
it was assumed that the drift mobility of the electrons was about 50 cm? 
volt-1 sec-1, the detailed analysis of the characteristics of the electrical 
circuit showed that their transit time would have been 5X 10~° sec. 
The fact that no latent image was ever formed on the rear surfaces of 
the pure crystals showed that, if electrons were actually liberated, their 
mean lifetime was considerably less than 5x10-®sec. This was in 
striking contrast to the lifetimes of 10sec which had been observed 
in air annealed crystals of silver chloride by Haynes and Shockley (1951). 
As already stated, Hedges and Mitchell (1953) had experienced no 
difficulty in displacing the photoelectrons in silver bromide crystals 
containing silver oxide in solid solution by the application of electric fields ; 
the lifetimes of the photoelectrons in these crystals were probably of 
the same order as those measured by Haynes and Shockley (1951) and 
Allemand and Rossel (1954). 

The crystals of silver halides which had previously been used for 
observations of photoconductive phenomena had always been grown 
in air and annealed in air at a temperature of about 400°c (Hecht 1932, 
van Heerden 1945, 1950, Haynes 1948, Allemand and Rossel 1954, 
Brown 1955). It was believed that air annealing (1) removed internal 
strains introduced during crystal growth and preparation of specimens 
and (2) removed the photochemical products responsible for the dis- 
coloration produced by exposure to light and thus restored the crystal 
to its original state (see, for example, Hecht 1932, van Heerden 1950, 
Allemand and Rossel 1954). When the reasons for the increased range 
of the electrons in the crystals which were grown and annealed in air 
were discussed, emphasis was placed on the removal of internal strains by 
annealing (Haynes 1948, Brown 1955) and the possibility that the crystals 
might be sensitized by chemical changes occurring during annealing was 
not considered. A crystal of pure silver halide which has been exposed to 
light to liberate halogen and to form particles of colloidal silver can, 
however, never be restored to its original state by annealing in air; the 
colloidal silver is oxidized to silver oxide which passes into solid solution 
in the silver halide. An exposed crystal of silver halide can be restored 
to its original state by annealing in air only if it was already a solid 
solution of silver oxide in the silver halide before exposure. Annealing 
in air does, however, remove positively charged groups of silver atoms 
which furnish traps for conduction electrons (Mitchell 1957) and con- 
sequently decrease their lifetimes and limit their ranges (Lehfeldt 1935). 
It is now clear that crystals which were in fact dilute solid solutions of 
of silver oxide in the silver halides have been used for all the experiments 
on photoconductivity and on crystal counters which have been recorded in 
the literature. Brown and Waifan (1957) have recently observed that 
the ranges of the electrons which are liberated by y-rays are greater in 
silver chloride crystals grown in air, and annealed in air at a temperature 
of 400° than in crystals grown in vacuum and annealed in vacuum, 
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in an inert atmosphere or in chlorine. Brown (1957) has also found 
that crystals of silver chloride, grown and annealed in air, have a greater 
photoconductive response when illuminated with short duration light 
flashes at room temperature. 

We are of the opinion that dissolved silver oxide plays an essential 
role in the photoconductive process which has not previously been 
recognized and we shall propose an interpretation of the high photo- 
conductive sensitivity of the air sensitized crystals. At low temperatures, 
crystals of silver halides containing silver oxide in solid solution have 
a higher ionic conductivity than pure crystals. This suggests that at 
least some of the silver oxide is dissolved in the form of doubly charged 
oxide ions, which occupy halide ion lattice sites, together with inter- 
stitial silver ions which compensate the additional negative charges ; 
the increased ionic conductivity is due to the drift of these interstitial 
ions. The primary process in electronic conductivity is the liberation 
of electrons and positive holes by photons or by charged particles. We 
assume that the lifetime of the positive hole is less than that of the 
electron. In silver halide crystals of high purity and perfection, the 
positive holes are trapped at the external surfaces and at dislocations. 
The resulting positively charged centre has, on account of Coulomb 
interaction, a large cross section for the capture of an electron so that 
the lifetime of a conduction electron is correspondingly short. Unless 
the positively charged centre dissociates with the loss of a silver ion into 
an interstitial position before recombination occurs (see Mitchell 1957), 
the external surfaces and the dislocations furnish recombination 
centres for positive holes and electrons. This is the reason for the small 
initial sensitivity of pure crystals of silver halides for the production 
of photocurrents and current pulses on irradiation. 

The air sensitized crystals have different properties. Here, the doubly 
negatively charged oxide ion has a large cross section for the capture 
of a positive hole, on account of the Coulomb interaction, whereas the 
resulting singly charged oxide ion has a very small cross section for the 
subsequent capture of a conduction electron. The interstitial silver 
ion which then compensates the negative charge of the conduction 
electron does not provide a trap for it except, possibly, at extremely 
low temperatures. After the positive holes have been trapped by doubly 
charged oxide ions they have to be released by thermal activation or 
optical excitation from the singly charged oxide ions before they can 
be trapped at the surface recombination centres. The result of dissolving 
molecules of silver oxide in crystals of silver halides is therefore to reduce 
the rate at which positive holes are captured by the surface traps, which 
act primarily as the recombination centres. The consequence of this is 
that the lifetime of the electrons in the conduction band is increased in 
the air sensitized crystals compared with crystals of the pure silver halide. 
This accounts for their increased photoconductive and photographic 
sensitivity. The correspondingly increased electron ranges provide an 
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explanation of why experiments of the type introduced by Haynes and 
Shockley (1948, 1951) may be successfully carried out with air sensitized 
crystals but fail with crystals of pure silver halides (Hvans 1955). 

This model for air sensitized crystals of silver halides which introduces 
two correlated energy levels in the forbidden band—an upper donor 
level associated with doubly charged oxide ions, and a lower acceptor 
level associated with singly charged oxide ions—should be compared 
with the double acceptor model for n-type germanium sensitized with 
iron, nickel, cobalt or manganese which has been developed by Newman, 
Tyler and Woodbury (Newman and Tyler 1954, Tyler, Newman and 
Woodbury 1955 a, b, Tyler and Newman 1955, Newman, Woodbury and 
Tyler 1956). The close parallelism will be immediately evident. Both 
models provide examples of the basic mechanism for sensitized photo- 
conductivity which has been outlined by Rose (1956). A further parallel 
may exist between the scattering of the conduction electrons by the 
doubly charged oxide ions in the sensitized silver halides (see the 
experimental results of Allemand and Rossel 1954, and of Brown 1955) 
and by the doubly charged manganese centres in germanium (Tyler 
and Woodbury 1956) as deduced from measurements of mobility at 
low temperatures. This interpretation of the origin of impurity scattering 
in air sensitized silver halide crystals at low temperatures is consistent 
with the estimated densities of impurity centres (Allemand and Rossel 
1954) and it seems likely that charged impurity scattering will always 
accompany air sensitization. 

The study of the temperature dependence of the electronic component 
of the electrical conductivity of solid solutions of silver oxide in silver 
bromide which are in contact with bromine vapour should yield informa- 
tion on the height above the full band of the lower acceptor level. The 
oxide ions will then be in the singly charged state and the crystal should 
behave as a p-type semiconductor because of the transference of electrons 
from the full band to the acceptor levels by thermal activation. From 
measurements at low temperatures of the wavelength dependence of the 
photoconductive response of air sensitized crystals of silver halides, it 
should be possible to determine the depth below the conduction band of 
the donor levels associated with the doubly charged oxide ions. 

The secondary photocurrents observed in air sensitized crystals at low 
temperatures by Lehfeldt (1935) may also be interpreted with this model. 
Under illumination, the upper levels in the forbidden band are emptied 
when the positive holes are trapped by doubly charged oxide ions. The 
lower levels associated with singly charged oxide ions which then appear 
can neither accept electrons from the full band at low temperatures nor, 
on account of the small cross section, capture electrons from the conduction 
band. In the steady state, a concentration of mobile conduction electrons 
with a compensating positive charge of immobile interstitial silver ions is 


therefore built up and comparatively large electronic dark currents can 
flow through the crystals. 
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This model thus resolves most of the difficulties which have arisen in 
the discussion of photoconductivity in the silver halides and of their 
behaviour as crystal counters. Their high photoconductive sensitivity, 
which has always been contrasted with the absence of photoconductivity 
in non-sensitized alkali halides, is due partly to sensitization and partly 
to intrinsic properties. It is clearly important to establish the intrinsic 
photoconductive behaviour of the silver halides over the whole range of 
temperatures ; measurements made by Braun and Mitchell (1957) 
with crystals of high purity silver bromide at temperatures between 
— 183° and —70°c will be published shortly. 

The presence of positive hole traps with the properties which have been 
attributed here to dissolved molecules of silver oxide appears to be a 
necessary condition for high photochemical sensitivity in silver halide 
crystals (Mitchell 1957, Mitchell and Mott 1957). 
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The Thermal Instability of a Rotating Fluid Sphere 
Heated Within: II} 
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SUMMARY 
It is shown that the characteristic value problem specifying the state 
of marginal stability of a rotating incompressible fluid sphere heated 
within, can be reformulated in terms of a variational principle. 


§ 1. INTRODUCTION 


In an earlier paper (Chandrasekhar 1957, this paper will be referred to 
hereafter as I) the problem of the thermal instability of an incompressible 
rotating fluid sphere has been considered in case the motions and the 
associated perturbations have symmetry about the axis of rotation ; and 
the underlying characteristic value problem was solved for a slightly 
modified set of boundary conditions. In this paper, it will be shown that 
the mathematical problem (as in other cases of thermal instability which 
have been investigated) can be reformulated in terms of a variational 
principle. This principle will be used in a later paper to solve the basic 
problem more rigorously than was accomplished in I. 


§ 2. THE VARIATIONAL PRINCIPLE 


The basic equations governing the problem are (I, eqns. (6), (16), (17) 
and (22)) : 


0 Aste) 
220V yfa 20 
42U — — —=+(—— — — é 
ee ey, “(5 =e Mee tt) 
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22 0U 
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where the various symbols have the same meaning asin I. By measuring 
distance in units of the radius, R, of the configuration and making the 


transformations 
2 
0— (= e)6 and v—-(=R) Vv, et ae A} 
Vv 
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eqns. (1)-(3) can be reduced to the more convenient forms 
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where 7’ and C are the Taylor and Rayleigh numbers, respectively, 
(cf. I, eqn. (38)). 

Multiply eqn. (5) by @ and integrate over the volume of the unit three- 
dimensional sphere. Both sides of the equation can be reduced by 
integration by parts and we find : 


i | | grad 6 [°r? dr du= { es eA =n) ur dee ee) 


the integrated parts vanishing on account of the boundary condition on 6 
(I, § 3.1). It will be observed that on the right-hand side of eqn. (8) 
the integration is over the unit five-dimensional sphere. 

Making use of eqn. (6) we can rewrite eqn. (8) in the form 


Of [erada prt dr du—[ [U(420—2 ) w? de de oe) 


Considering first the term in U on the right-hand side of this equation, 
we can reduce it by making use of eqn. (7) and a sequence of integrations 
by parts. Thus, 


— [|S oder de 
Oz 
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the integrated parts do not make any contributions in view of the 
conditions that U and either V or dV/dr vanish on the boundary (cf. I, 
eqns. (34) and (35)). 
Next, letting 
A=4,U, RP gee tt? ere tae LL) 
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we can transform the remaining integral in eqn. (9) as follows : 


| | UA,Xa? dew dz 


=| ‘U {=p = G Ge) + 3 a =| =|} dr dys 
=—|{ {ro a ua ee Ue Nae eh 
=-{" (F =) (X},10—w? ) du+| | XA,Ur(1—p2) drdu. (12) 


On a rigid boundary dU/dér=0 (I, eqn. (30)) and the ‘ surface ’ integral 
in the last line of eqn. (12) a cece On the other hand, on a free 


boundary (I, eqn. 32)) 
CUT 2 OU, 


ore | r Or oo 
Consequently, 
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Finally, combining eqns. (9), (10) and (15), we obtain, 


alee +1 T\2 
c| | | grad 6 Pr? dr du=—2|_ (S) (1—p2) du 
(ae) ail or rei 


ss 
+) J {4s i +1| (5) 3 F(z) |ptane dr du. (16) 


This is the required formula which, it can be readily shown, provides the 
basis for a variational procedure for determining C. 

One general result} which can be proved from eqn. (16) is that CO is a 
monotonic increasing function of 7. 
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CORRESPONDENCE 


On Spiral Etch Pits in Germanium and Silicon 


By 8. G. ELuis 
RCA Laboratories, Princeton, New Jersey 


[Received July 16, 1957] 


Ty an earlier communication (Ellis 1955) I have reported the observation 
of etch pits with a spiral terrace in germanium crystals. They were there 
ascribed to the nucleation of the etching by screw dislocations with Burgers 
vectors of the order 10-° cm. Although screw dislocations of this strength 
were not inconceivable in the dendritic crystals studied, the observation 
of similar pits in germanium and silicon crystals grown slowly from the 
melt casts doubt on the above interpretation. Etch pits of this sort 
have now been observed in a number of cases (Christian and Jensen 1956, 
Bardsley and Straughan 1956, Vogel and Lovell 1956), and Lang (1957) 
has suggested that they are formed when the rate of etching parallel to 
the surface of the crystal is asymmetric about the point in the surface 
at which the attack is nucleated. 

A simpler explanation is suggested by some recent observations of 
Bontinck and Amelinckx (1957). By decorating the dislocations in 
synthetic fluorite crystals they observe that some of the dislocations form 
helices when a dislocation of mixed character (i.e. a composite dislocation) 
climbs as vacancies condense on it. If such dislocations exist in ger- 
manium and silicon and if the axis of the helix is nearly normal to the 
surface being etched, then the terraces in the pit will spread from a point 
which has moved in a circle as the etching proceeds. Figure 2 of Lang’s 
paper then shows how with the limited resolution of the light microscope 
a spiral terrace appears. This explanation differs from Lang’s in supposing 
that the lateral rate of etching is constant and that the point of nucleation 
moves in a circle. 

If this explanation is correct the lateral distance between terraces 
is a measure of the diameter of the helix, and the step height between 
terraces is a measure of its pitch. The pitch, as Bontinck and Amelinckx 
point out, can be much larger than the Burgers vector of the dislocation. 

While these ideas are consistent with the observations of etched crystals 
such observations cannot give complete confirmation. This would be 
better done by decoration methods. If the idea is correct it provides 
further information about the character of the original dislocations, and 
the extent of vacancy condensation at them. 

In conclusion it may be added that spiral etch pits have also been seen 
on the (100) and (110) surfaces of germanium crystals. 
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REVIEWS OF BOOKS 


Progress in Cosmic Ray Physics, Vol. III. Edited by J. G. Witson  (North- 
Holland, Amsterdam.) [Pp. xii+420.] Price 76s. 


Tuts third report on recent important advances in cosmic ray physics contains 
valuable articles by K. Greisen on Extensive Air Showers; H. 8. Bridge on 
Charged K-mesons and Hyperons; R. W. Thompson on Neutral K-mesons 
and Hyperons, and G. Puppi on the Energy Balance of Cosmic Radiation. 
This volume, published late in 1956, is no more up to date than Spring 1955 ; 
strange particle physics grows and changes rapidly, so this delay is unfortunate. 
It is intended that future volumes will cover elementary particle theory as 
well as cosmic ray physics. J.-H: 


Quantum Mechanics. By H. A. Kramers, English translation by D. TER Haar. 
(Amsterdam : North-Holland), [Pp. xvi+496]. Price 90s. 


Dr. TER HAAR has made a very satisfactory translation of Kramers’ well known 
monograph on quantum mechanics. The original was completed in 1937 ; the 
scope of the work does not go much beyond the quantum theory of atomic 
systems. spas fe 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. ] 


S. AMELINCKX etal. Phil. Mag. Ser. 8, Vol. 2, Pl. 48. 


Etch.pits on (111) cleavage faces of Cal’,. («) Observe the bunching of steps 
in the [211] direction. (4) Etch-pits having the sector of bunched steps 
in different [112] directions. 


Fig. 6 


Pik 


Etch-pit due to a sequence of prismatic loops. Notice the radial discontinuity 
in a [112] direction. 


